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MASS SPECTROMETER CALCULATIONS 
on the 
IBM 602-A CALCULATING PUNCH 



W. H. King, Jr. and William Priestley, Jr. 
Esso Laboratories, Standard Oil Development Company 

INTRODUCTION 

The following instructions describe in detail the control panels for the IBM Type 602-A 
Calculating Punch which are used to perform a matrix multiplication and subsequent 
normalization steps. The problem, which has been discussed in detail in another 
paper 1 , is discussed here in brief for reference purposes. 

THE PROBLEM 

1. To multiply a 20th order matrix by a vector. 

2. To calculate from the answers obtained in (1) above the following percentages: 

a. Total normalized mole per cent. 

b. Air free mole per cent. 

c. Fixed gas free mole per cent. 

d. Fixed gas free weight per cent. 

e. Air free specific gravity. 

f. Fixed gas free specific gravity. 

The inverted matrix shown in Figure 1 will be multiplied by the peak heights 
(b , b , b , etc.) to obtain 19 unnormalized mole per cent answers (X x , X 2 , X 3 , etc.) 
and a summation check. Since a summation check is calculated along with the 19 
answers one has effectively a 20 x 20 matrix multiplication problem. To obtain any 
given X value, for exampleX^ one multiplies the peak heights b x , b 2 , b 3 , etc. by their 
corresponding matrix values B, ;B ;B ;etc. A summation of these products yields 
X . Similarly X 19 is obtained by multiplying the peak heights by their corresponding 
matrix values for the 19th row B 19 l ; B 19 2 ; B 19 ; etc. In the twentieth row the B n 
values are the sum of the columns,' i. e. , B n 5 is the sum of all B's in column 5 
(B + B + B + etc. ). The twentieth row is therefore the summation check and 

'1,5 2,5 35' 

the value obtained for Z X n should equal the sum of X A + X 2 + X 3 + X 4 +. . . + X lg if 
the matrix has been multiplied correctly. Since 2 X n is the sum of 19 multiplications 
and the value from totaling X 1 + X 2 + X 3 etc. is the sum of 361 multiplications, a com- 
parison of the two figures is an adequate check on the machine's computation. 

* See authors' article, "Spectrometic Analysis Employing Punch Card Calculators", 
in Analytical Chemistry, Vol. 23, No. 10 (October, 1951), pp. 1418-1420. 



It was decided to make the matrix multiplication in blocks of 8, as indicated in Figure 
1, in order to take advantage of the time saved by making two simultaneous multiplications. 
This calculation method also utilizes the standard IBM 602-A Calculating Punch at maxi- 
mum capacity. (Additional capacity is available at increased cost.) Two rows of the 
matrix are calculated simultaneously, resulting in two answers for each five matrix 
cards. The matrix has been so arranged that the 17, 18, and 19th answers are fixed 
gases N 2 -CO, C0 2 and air. 

The normalization procedure consists of determining the normalization factors and 
making the division steps in order to provide the required answers. 

These factors and steps are symbolized below: 



(Total sample mole % normalization factor) 
1. X 1 + X 2+ X 3+ ... + X l9 = S M% 
(Normalization factor for the mole % air free) 



2. X 1 + X 2 + X 3+ ... + X l8 =SM% 



i 

s 

18 

(Normalization factor for the mole % fixed gas free sample) 

3. X 1 + X 2+ X 3+ ... + X 16 =SM% 

(Normalization factor for the weight % of fixed gas free samples and unnormalized 
specific gravity for this sample) 

4. (X 1 )(M.W.F. 1 ) + (X 2 )(M.W.F. 2 ) + ...+ (X 16 )(M.W.F. 16 )= | Pr 
(Unnormalized specific gravity for the air free sample) 

5. (XJ (M. W.F.J + (X 2 ) (M.W.F. 2 )+... + (X l8 )(M.W.F. ) = Z Pr 



18 



Molecular Weight Factor: 

6. (M. W. F.) = molecular weight of the component 

molecular weight of air 

Total mole per cent answer: 

100 X , m 

7. M% t = = unnormalized mole % x 10Q 

2 - M % . normalization factor 

19 

Air free mole per cent answer: 

100X 

8. M% af = 



S M% 

19 



Fixed gas free mole per cent answer: 

100X 

9- M% fgf = r 



S M % 

16 



Fixed gas free weight per cent answer: 
10. Wt.%. f = (X) (M.W.F.) (100) 







S Pr 

16 


Aii 


• free 


specific gravity: 


11. 


S.G. 


, f = (S Pr) (100) 

18 




S M% 

18 


Fixed gas free specific gravity: 


12. 


S.G, 


r.f = (S Pr) (100) 



S M% 

16 

For those not familiar with the actual equipment, a brief explanation of the cal- 
culating punch will be of interest. The 602-A has three basic operations; reading in- 
formation from a card, performing calculations, and punching results. Information 
read from a card may be entered into several storage units designated as 1L, 1R, 2L, 2R, 

3L, 3R, 4L, 6L, 6R, 7L, 7R. A number existing in a storage unit will remain there until 
another number is read into that storage unit. Numbers may be read out of any stor- 
age unit many times for calculation purposes without destroying the number. Storage 
unit 1R is unique since divisors or multipliers must be entered here (from other stor- 
age units or from reading of the card) just prior to the division or multiplication. The 
standard 602-A also has 6 counter units where products, quotients, or sums may be 
developed and accumulated. There, counters may also be used for storage. 



MATRIX CALCULATION CONTROL PANEL #1 
(Figures 2 and 3) 

The problem at hand requires the storage of 12 numbers, 8 matrix values and 4 
multipliers. The planning chart F 182 has been filled out for the first block of 8 
multiplications. The 602-A will perform two simultaneous multiplications by a com- 
mon multiplier resulting in two ten digit products. It is possible to perform more 
simultaneous multiplications if the number of required digits is less. It is advantageous 
to do this in order to save time. In this control panel setup, thesumofproductsoftheodd 
rows will be accumulated in counters (2, 3) where they are developed. Simultaneously, 
counters (5, 6) will develop and accumulate the products of the even rows. Thus, when 
the 20th matrix value has been calculated for the first two rows, counters (2, 3) will 
contain 



M B 1?1 ) + »,<%> + b s<B lj3 ) +-.- + hu,<B lfl9 ) 



and counters (5,6) will contain 

MB 2 , J + b 2 (B 2> 2 ) + b 3 (B 23 ) + . . .+ b l9 (B 219 ). 

These two sums will be the first two answers and will be punched on trailer cards fol- 
lowing the five detail or matrix cards. 

Read Cycle 

On the read cycle, the multipliers b 1 , b 2 , b and b 4 are read from the matrix 
card into storage units 1R, 2L, 2R, 3L, while the multiplicands B , B , B , 
etc. are read into storage units 3R, 4L, 4R, 6L, 6R, 7L, 7R and counter (1) . The 
numbers under the matrix symbol in Figure 2 indicate their column position in the card. 
For example B 2 1 is a 4 digit number in columns 62, 63, 64, 65 while column 61 is a 

or 5 to indicate the sign of the 4 digit number. 50064 means - 64, while 00064 means 
+ 64. After the read cycle, the calculator proceeds to the program steps. 

Program Step 1 (For Sign Control) 

Since the first multiplication will be (\ ) (B x 1 ) and (\ ) (B 2 1 ), the multiplier was 
entered in 1R on the read cycle. The sign control for the multiplication is usually ac- 
complished on the same program step as the multiplication by using a or 9 in a speci- 
fied digit position of the multiplicand. This would require, for 8 multiplications, the use 
of 8 pilot selectors (sign controllers). The standard 602-A has 7 pilot selectors and con- 
sequently a different system must be employed. By using an extra program step per 
multiplication, two pilot selectors may be used four times to accomplish sign control. 

In this program a or 5 is read out (R. O.) of the first digit position of storage 
unit 3R, and the board will be wired to permit a 5 to energize pilot selector number 

1 to indicate that the sign of B is negative. If a exists in this digit position, 
pilot selector 1 will not be energized (indicating B is positive) and the succeeding 
multiplication will proceed positively. A similar arrangement is used for the sign con- 
trol of B , , via pilot selector number 2. These two pilot selectors will be used over 

2,1 

and over as will become apparent from the succeeding discussion. Pilot selector num- 
ber 1 (P.S. #1) will always control the sign of all odd row matrix values, and (P.S.#2) 
will control the sign of all even row matrix values. No sign control is required for the 
vector or multiplier since in this problem the value is always positive. These pilot 
selectors must be reset to clear out the previous sign control. This is accomplished 
by dropping out (resetting) the pilot selectors with the digit selector before a 5 from the 
storage unit can energize it. In the operation of the machine, digits are read in a time 
sequence, 9, 8, 7, 6, 5, 4, 3, 2, 1. (No impulses are emitted for a zero.) When the 
digit selector is impulsed, it immediately emits a series of impulses in this order. The 
first impulse (9) is wired to drop out the selectors. This occurs before the impulse for 
(5) from the storage unit can arrive at the pilot selector for sign control. 

Program Step 2 

Multiplication is accomplished by impulsing the multiplier. In this operation the 
value in storage unit 3R is read into counter (2, 3) through the sign control pilot selector 



for as many times as there are digits in the multiplier 1R. For example, if (1R) is 123 
and (3R) is 00366, the value 366 will enter counter (2, 3) 123 times. In reality 36600 
will enter once, 3660 twice, and 366 three times. The accumulation therefore will be 
40626. If 3R were 50666, indicating minus 366, 36600 would be deducted from counter 
(2, 3) once, 3660 deducted twice, and 366 three times. This would leave 9999959374 in 
the counter which is the complement of 40626; the 9 in the extreme left digit position 
serves to indicate it as a complement, not a positive value. 

Program Step 3 

This step is similar to program step 1 (P x ) except that the multiplier is changed 
to b by reading b out of storage unit 2L into the multiplier storage unit 1R. The sign 
control of B and B is accomplished by reading the fifth digit positions of (4L) 

1)2 « j 

and (7L) to their respective pilot selectors. 

Program Step 4 

Multiplication here is again similar to P 2 . The counters will now contain the sum 
of the products of the first two sets of multiplications. Program steps 5, 6, 7 and 8 are 
self-explanatory since they are a repetition of the early programs which have been ex- 
plained in detail. At this point, the data furnished on detail or matrix card # 1 have 
been calculated. 
Program Step 9 

Since counter (1) was used as a storage unit it must be cleared before another num- 
ber is entered, hence, counter (1) is reset on this program. The sign control pilot 
selector is also reset by impulsing the digit selector on this program. The next opera- 
tion of the machine will be reading a new matrix card, consequently, the read hub is im- 
pulsed. Programs & 1 through *# 9 will be repeated for each of 5 matrix cards. 

After the 5th detail or matrix cards have been computed, counters (2, 3) will contain 
the odd row answer while counters (5, 6) will contain the even row answer. 

Program Step 10 

When a trailer or answer card enters the machine with an "X" punch in column 14, 
all preceeding programs are skipped and the programming will start here. 

The function of X14 cards will be to punch the answer of odd rows in counters (2, 3) 
in this card. Counter (2, 3) is read into punching storage unit 6 and the answer will be 
punched from this storage unit in the card. 

The next step is to read a new card; consequently read is impulsed. 

Program Step 11 

When a card enters the machine with an "X" punch in column 15 all programs pre- 
ceeding program 11 will be skipped and programming will start here. The function of 



This is in reality a simplification of how the IBM Type 604 Electronic Calculating 
Punch multiplies. The actual multiplication procedure for the Type 602-A is much 
more involved. 



XI 5 cards will be to punch the answer for even rows in this card. 

Summary of Operation Thus Far 

Five detail (or matrix) cards have entered the machine and have been calculated. The 
results of the calculations have been punched on two trailer cards. The answer to the 
first row multiplications is in the X 14 trailer card and that of the second row in the X 15 trailer 
card. 

This process is iterated for row 3 and 4 etc. until the entire matrix has been com- 
puted. After the matrix has been computed there will be 70 cards to sort. Fifty of these 
are the matrix or detail cards, while the remaining 20 are the X14, X15 trailer cards 
containing the 19 unnormalized mole per cent answers and the check row answer. The 
next control panel will process these 19 answer cards and produce three more cards containing 
the normalization factors. 

The system used here is not limited to a 20 x 20 matrix. Any order matrix may be 
used by extension of the cards entering the machine prior to the trailer cards and no 
control panel changes are required. 

Figure 3 shows the wiring diagram for the matrix calculation. This drawing is the 
actual wiring involved in the control panel. 



MOLECULAR WEIGHT FACTOR MULTIPLICATION 
AND SUMMATION CONTROL PANEL #2 

Figures 4 and 5 show the planning chart for the molecular weight factor multiplication 
and the summation of certain answers and products. 

Answer cards (X14 and X15) were punched on control panel#l when the matrix was 
calculated. In these cards the corresponding molecular weight factor must be punched 
by a gang punch or hand punch prior to this step. 

The function of this control panel is: 

(a) To multiply the 19 unnormalized mole % answers by the molecular weight 
factor. (M.W. F.) 

(b) To determine the normalization factors: 

1 • C S M% ) the total sample mole % normalization factor which is the sum 
of the 19 unnormalized mole per cent answers. 

2. ( S M% ) the air free sample mole % normalization factor which is the 

18 

sum of 18 unnormalized mole % answers (omitting the air card). 

3. ( S Pr ) the fixed gas free sample mole % answer which is the sum of 16 

16 

unnormalized mole % answers (omitting air, nitrogen-carbon monoxide, 
and carbon dioxide) . 

i 

4. (S Pr) the fixed gas free weight % normalization factor which is the sum 

of the products of the mole % times the molecular weight factor for 16 



10 



cards. This factor when divided by 2 M% yields the specific gravity for 

16 

the fixed gas free sample. 
5. ( 2 Pr ) the sum of the mole % times the molecular weight factor for 18 

18 1 -_, 

cards. This factor when divided by 2 M % yields the specific gravity 

18 

for the air free sample, 
(c) To punch these results in the appropriate cards. 

Read Cycle 

The molecular weight factor (M.W. F.) in card columns 25 through 29 is read into the 
multiplier storage unit 1R; the unnormalized mole % answer in columns 17 through 21 is 
read to both storage unit 2R and counter (1). This will accumulate the sum of the un- 
normalized mole % answers. 

Program Step 1 

The unnormalized mole % answer (M %) is multiplied by (M.W. F.) by reading 2R 
into counters (5, 6), impulsing multiply, and reading out 1R. 

Program Step 2 

Because the products obtained in program step 1 are to be totaled, this program is 
used to transfer the product to counters (2, 3). At the same time the product is read into 
punching storage. 

Program Step 3 

This program punches the product in this card and sets the machine to read the next 
card. Programs 1 through 3 are repeated for each of the 19 unnormalized mole % cards. 

Program Step 4 

When an answer card containing an "X" in column 22 enters the machine, programs 
1 and 2 are performed but program 3 is skipped. 

At this time 16 cards have been computed and therefore counter (1) contains the sum 
of the uncorrected mole % answers ( 2 M % ) thus far. Likewise, counters (5, 6) 
contain the sum of the products (2 Pr). 

16 
1 1 

2 M % is stored in 3L and 2 Pr is stored in 3R at this time. These results will 
be punched later on trailer cards. 

Program Step 5 

P is identical to P and the product is punched in this card (X 22) and "read" is 
impulsed to proceed with the next card. 

Program Steps 6, 7, 8 

These steps are similar to the above except that 2 Pr are entered in 4L and 4R 

1^1 18 

while 2 M % and 2 Pr are left in the counters. At this time 19 answer cards have 

19 19 

been processed and the normalization factors have been placed in storage. The next few 
program steps are used to punch the normalization factors on three trailer cards. 
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Program Step 9 

When a trailer card with an "X" in column 52 enters the machine, programs 1 
through 8 are skipped. In this program s M % is read out of 3L and read into 7R for 

* 16 

punching. Similarly, 2 M%xM.W.F. enters 6R. 

IB 

Program Step 10 

Both punching storage units are impulsed to punch the summations in the X52 card. 
Read is impulsed to start the next card. The X52 card now contains the normalizing 
factors ( SM%) ( 2 Pr) for the final computation on the next control panel. 

*6 16 

Programs 11, 12, 13, and 14 

These programs are similar to programs 9 and 10 with the exception that they are 
initiated by X55 and X60 cards. 

Summary of Cards Thus Far 

Nineteen answer cards have been processed and three trailer cards have been pro- 
duced containing the partial and total summations. 

Each answer card now contains 

A. Unnormalized answer 

B. Molecular weight factor 

C. A times B 

X52 card contains 

2 M % and 2 Pr. These are the normalizing factors for obtaining an 

16 16 

analysis without air, nitrogen-carbon monoxide and carbon dioxide. 

X55 card contains 

1 i 

2 M % and 2 Pr. These are the air free normalization factors. 

18 18 

X60card contains 

1 i 

2 M % and 2 Pr which are the total sample normalization factors. 

19 19 

Figure 5 is the actual control panel wiring for the above calculations . Normally the 
602-A Calculating Punch has 12 programs, however, 14 were required for this opera- 
tion. The two extra program steps were obtained on this problem by wiring the program- 
exit hubs of P x and P through co-selectors #■ 5 and & 6. When the co-selectors are 
not energized, (the normal position) P A andP 2 are performed. When a card with "X" 
punched in column 60 enters the machine the co-selectors are energized by pilot selector 
#" 6. P 1 and P 2 when wired from the transfer side of the co-selectors actually perform 
P. andP,^ instead of P , and P„ . In other words, P, andP„ exit hubs are used 

13 14 12 '12 

twice, once to perform P x and P 2 and once for P andP 14 - 
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THE NORMALIZATION CONTROL PANELS 3 
(Figures 6 and 7) 

Read Cycle 

The read cycle has been split into four parts in order to enter the divisors from 
X52, X55, and X60 cards into storage. 

The X52 card picks up pilot selector#l through the control brushes and reads 
S Pr into storage unit 6L and 2 M % into storage unit 6R. Similarly, the X5 5 
card reads E Pr and 2 M % into 2L and 2R, while the X60 card reads S Pr 

1 18 18 19 

and S M % into 3L and 3R. 

19 

The next card to be read will be an X14 or an XI 5 answer card. These cards will 
pick up pilot selector # 4 to control the reading of the unnormalized mole per cent 
answer (M %) into storage unit 4R and the product M % X M. W. F. into storage unit 4L. 

Program Step 1 (D i vision Set Up) 

The divisor 2 M % in storage unit 6R is read into the division storage unit 1R 
and the dividend M*^, in storage unit 4R is read into the dividend counters (1, 2, 3). 

Program Step 2 (Division) 

In order to perform division the divisor is subtracted from the dividend in a se- 
quence until the dividend is reduced to zero. The number of times required to do this 
will be the quotient which will be developed in the counters (5, 6). 

1R is read into (1, 2, 3) negatively, the quotient is read into (5, 6) and the divide 
hub is impulsed. 

Program Step 3 (Punch Normalized Answer) 

Counters (5, 6) now contain the normalized answer for M % fixed gas free. This 
value is punched. 

Counters (5, 6) are read out to punching storage 7R and punched in the card. At the 
same time (5, 6) are reset to 5 in order to round off the answer. At this time the card 
is held in the machine for the next division. 

Program Steps 4, 5, 6 

These steps are similar to P x , P 2 , P 3 and the normalized M % air free answer is 
obtained because the divisor this time was S M % . 

18 

Program Steps 7, 8, 9 

These steps are again similar to P x , P 2 and P 3 and the total M % normalized 
answer is obtained by using 2 M % for the divisor. 

19 

Program Steps 10, 11, 12 

These steps calculate the normalized weight per cent fixed gas free answer by chang- 
ing the dividend to M % X M. W. F. and the divisor to 2 Pr. 

At the completion of program 12 the read hub is impulsed through the normal side of pilot 
selector#6, and the machine reads the next card. 
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Program Step 13 

The next six program steps will be initiated only for answer cards containing an "X" 
in column 28. The X28 card will pick up pilot selector *t 6 and after program step 12 has 
been completed, program 1 exit hubs will again become active because the read hub im- 
pulse from program 12 will be broken by pilot selector ^ 6. At this time, co-selectors 
1, 5 and 8 will be energized and program 1 will function as program 13. 

In this program S Pr in storage unit 6L is read into the dividend counters (1, 2, 3) 
and S M % in storage unit 6R is read into the divisor storage unit 1R. 

Program Steps 14 and 15 

These steps are identical to P 2 and P 3 except that specific gravity fixed gas free is 
computed and punched. 

Program Steps 16, 17, 18 

These programs are similar to programs 13, 14, and 15 and the air free specific 
gravity is calculated and punched. 

Summary of Cards Processed Thus Far 

1. The matrix calculation produced 19 answer cards. 

2. The summation and multiplication calculation punched the M % x M. W. F. in the 
19 answer cards and produced three cards containing the six normalization 
factors. 

3. The normalization calculation punched four normalized answers in each of the 
19 answer cards and in the 19th answer card and the two specific gravities were 
also punched. 

Figure 7 shows the control panel for the normalization control panel. In 
general, the wiring is straightforward and may be duplicated from the planning chart 
by persons familiar with the equipment. Programs 2, 5, 8, 11, 14 and 17 are in all 
respects identical and have been so wired. Programs 3, 6, 9, 12, 15 and 18 are also 
identical. On programs 1, 4, 7, 10, 13 and 16, storage unit 1 R and counters (1, 2, 3) 
are read in and have been so wired. One exit hub of program 1 has been wired through 
co-selector # 5 to have this exit hub read out of 4 R on program 1 and read out 6 L on 
program 13. Similarly, one exit hub of program 4 has been wired normally through 
co-selector # 5 to read out 4R. On program 16, this exit hub reads out 2 L. 

In order to prevent a back circuit with the exit wiring from storage unit 6 L, 
these wires were run through co-selector # 1. Co-selector #8 was used to stop the 
skip impulse from interfering with punching in column 61 when the specific gravities are 
punched in the X28 card. 
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COMPUTATIONS OF INVERSE MATRICES 
BY MEANS OF IBM MACHINES 



Jack Sherman 
The Texas Company Research Laboratories 



An efficient method for solving the sets of simultaneous linear equations which 
arise in mass spectrometer analyses is to utilize the inverse matrix of the calibration 
coefficients. Once the matrix has been obtained, its application to the solution of 
simultaneous linear equations merely consists of computing sums of binary products. 
However, the task of obtaining the inverse matrix is considerable inasmuch as the order 
of the matrix generally ranges from 13-20 in routine operations. 

Many methods have been and are being published for obtaining the inverse of a 

a 

matrix . These generally involve some variation of the method of systematic elimination, 
such as Doolittle's, and were designed with the idea that the computations would be 
carried out by means of a desk calculator. It turns out that these methods are not par- 
ticularly efficient for use with IBM machines. 

It is the purpose of this paper to describe an efficient method of obtaining an 
inverse matrix by means of punched cards which has been found to be relatively easy to 
carry out and which requires considerably less time than for any of the standard methods. 

The computation of the inverse of a matrix of pattern or sensitivity coef- 
ficients is begun by means of the previously obtained inverse (as will be seen in the follow- 
ing description of the method, the final inverse does not require that this be the starting 
point. However, it is convenient to do so.) This starting inverse is transformed into 
another inverse which results from replacing the first column in the original matrix by 
the first column corresponding to the new matrix. The equations for doing this are 
given in an abstract of the paper presented at the Twelfth Summer Meeting of the Institute 
of Mathematical Statistics at Boulder, Colorado. b This may be written as follows: 

a Cf.H. Hotelling, "Some New Methods in Matrix Calculations", Annals of Mathematical 
Statistics, 14, pp. 1-34 (1943). 

b Annals of Mathematical Statistics , Volume 20, page 621, (1949). 
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Denote the original inverse by [b] whose elements are b^, the new inverse 
by [b'] whose elements are b'. , the original matrix by [a] whose elements are 
a , and the new matrix by [a'] whose elements are a' . If [a'] differs from [a] 
only in the elements of the k th column, then 

b;. = b t - z f b; , i = l,2,...,k-l,k+l,...N (1) 

j = 1,2,... N 

Ki = b kj /z k , i = 1,2,... N (2) 

N 
z = y b, a' , i = 1,2, ...N (3) 

i £_ i r rfc ' ' ' 

r=l 

By letting k = 1,2, . . .N successively, the new inverse matrix is obtained from 
the old inverse matrix by successively replacing each column of the given matrix [ a ] by 
a column of the new matrix [a'] . The N successive matrices obtained in this process 
correspond to inverses of the old matrix [ a ] in which the first column has been replaced 
by the new calibration coefficients, the old inverse in which the first two columns have 
been replaced by new calibration coefficients, etc. 

Although the foregoing procedure of transforming the old inverse matrix into a 
new inverse appears to involve a greater absolute number of arithmetical steps than per- 
haps the Doolittle method, it can be very efficiently carried out on the IBM machines. The 
following series of steps which are utilized in making one transformation are listed below: 

1) N 2 cards are prepared by punching one element of the starting inverse into 
a card in columns 26-31, the order of the cards being elements of successive rows, 
proceeding from left to right. 

A negative value of an element is indicated by an Xin column 31. These cards 
contain anX in column 52 for control purposes. 

For identification purposes the values of i and j for each element and the 
number of the inverse are punched in columns 1-6, thus 050800 denotes b 5g of matrix 
0, the starting matrix. 

2) A second set of N 2 cards is prepared containing the values of the indices 
i and j in columns 1-4 together with the number 01 in columns 5 and 6 to denote 
inverse matrix number 1 (the [b '] matrix). 

3) The N values of a' are reproduced N times into N 2 cards in columns 
71-76 of matrix 0. 
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4) Behind each set of N cards of inverse matrix 0, a trailer card is inserted 
by means of the collator. 

5) The z i 's are computed utilizing equation 3. These z 's are punched in 
columns 60-68. An X in column 68 denotes a negative value of z. 

It should be noted that if the original elements of the first column of a are 
used rather than a new set of elements, the values of the z's would be 1, 0, 0, . . . . 
If the calibration coefficients a^. do not differ greatly from the original values a . , the 
z values will tend to be close to the 1, 0, 0, . . . , and consequently the order of 
magnitude of the z values is well controlled. This is the reason for starting with the 
old inverse. 

6) The reciprocal of z v is obtained by means of a desk calculator (or 
reciprocal control panel) and a new card is punched with 

z' = J- 

i — > 

the digits being placed in corresponding positions of the z * cards. 

7) The %\ value is gang-punched into a new set of N cards in columns 61-68, 
at the same time reproducing b } . values from card columns 26-31 of the first N cards 
of original matrix into card columns 16-21. 

8) The first N cards of matrix number 1 are interspersed with the N cards 
obtained in the previous step by means of the collator. 

9) The values of b' are computed in card columns 16-21 according to 
equation 2 which may be written as 

b ii = °+ z ; * lt , (4) 

the computations being arranged so that z[ and b are obtained from a card prepared 
as described in step number 7, and the corresponding values of b '. being punched on 
the following card, i.e., a card of the number 1 matrix. 

10) The values of b*. are reproduced into card columns 16-21 of the remaining 
N 2 -N cards of the [ b] matrix (matrix 0); at the same time the z values are gang- 
punched into columns 61-68. 

11) The corresponding N 2 -N b' cards (matrix number 1) are interspersed 
with the b cards by means of the collator. 

12) The b ' values are computed according to the equation 1. 

It should be noted that equations 1 and 4 are similar in form, the only differ- 
ence being that for equation 4 the product term is added whereas in equation 1 it is 
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subtracted. Hence, the same control panel can be used for both computations. 

13) The set of two (N 2 — N) cards are separated by means of the sorter into 
the b and b * cards. For a 13th order matrix, each of the N transformations can be 
carried out in approximately 25 to 30 minutes, and so the desired inverse can be 
obtained in approximately seven hours. 

The following points should be noted concerning the method of a computation 
described in the foregoing steps: 

A running check not described in the foregoing series of steps is introduced 
in order to detect any error before the computation is completed. This consists of 
inserting a column of the a' matrix into the cards so that when the z's are computed 
according to step 5, additional summations are carried out, -namely, £_/ b i r a r • 
These summations must yield the. values 0, 0, . . . 1, 0, . . . . 

After the foregoing steps are carried out N times to obtain the desired inverse, 
a complete check is carried out by multiplying the new inverse b ' with the correspond- 
ing given matrix a' to obtain the unit matrix. 

Difficulties associated with the position of the decimal point in division are 
avoided by the fact that the z' values are obtained separately by means of a desk 
calculator (or reciprocal control panel) . 

In addition to the advantage that is gained by making the. computation com- 
pletely routine, the method of inverting a matrix described in the previous steps provides 
a considerable degree of flexibility; thus, if only a few of the N components are re- 
calibrated, the corresponding new inverse can be obtained with considerably less comp- 
utation than would be required for the complete calibration. 

Although the series of steps described in this paper for computing an inverse 
matrix requires the use of the calculator, sorter, reproducer, and collator, the sequence 
is easily learned by the machine operator after one or two examples are carried out. 
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LINEAR EQUATIONS AND MATRIX INVERSION 

Eric V. Hankam 
Watson Scientific Computing Laboratory 



A. Mathematical Preliminaries 
1. Linear Equations: 

Consider a square matrix of order n (let n = 4 for simplicity) 

a i 2 a i3 a i4 b is 
a a a b 

22 23 24 25 

a 32 a 33 a 34 D 35 

a., a., a.. b„ 

42 43 44 45 

augmented by a column vector b, representing a set of four 
simultaneous linear equations. 

I. a) Perform the following transformation: Divide each element of 
the first row by its leading element (a u ), thus obtaining: 




1 


a !2 


a !3 


a i4 


b X5 




a n 


a n 


a il 


a 

ii 


a 21 


a 22- 






b 25 


a 41 








b 45 



(2) 



26 



b) Subract a times the first row from the respective elements 
of the succeeding rows (i = 2, 3, 4): 



a 22 " a 21 






(3) 



Note: In place of (2) each row may be divided by its leading elements: 



13 



21 



(2*) 



42 



The first row is then subtracted from all the other rows: 



12 



15 



1 


a 


11 






a 


ii 





a 22 


a !2 
a il 


• • 


b 2S 


- 


b 15 


a 21 


a 21 


a il 







. 






■ 







a 42 


' a i2 
a il 


. . . 


b 4 5 

a 

41 


- 


fa l 5 


a 41 


a il 



(3*) 



The reduced matrix of (3) or (3*) may be written in the form 



"24 



35 



(4) 



27 



II. Applying to (4) the same transformation as in step I. , (4) is reduced to 
a second order matrix: 



a 34 b 35 



' a 42 a 44 



III. Reducing (5) similarly, the following is obtained: 

Performing the transformation of step la) on (6), we obtain 
1 ¥- = 



(5) 



(6) 



(V) 



where x 4 = the fourth variable of the set of equations. 
This leaves the original matrix in a triangular form: 



12 



13 



23 



a ll 


a 24 


a 22 


a 34 



33 



15 



25 



K 



(8) 



4 5 _ 



= x. 



(8) 


may be written: 




/■ 


a 




[• 


1 




• 







\. 






12 



13 



23 



1 





24 



34 



15 



35 



b 45 =X < 



(8) 



IV. In order to reduce (8) to the unit matrix, the following transformations are 
performed on rows 1, 2, 3 of (8): Multiply an element of the fourth column 
by the elements of the last row and subtract the products from the respective 
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row yielding (9) (i. e. , a^ - a j4 a 4j ; i = 1, 2, 3 ; j =' 1, 2, . . . , 5). This trans- 
formation leaves columns 1, 2, 3 unchanged; 



1 


a 





1 









12 



13 



23 






b 15 _ a !4 X 4 







"b 25 _ ^24 X 4 







b 35 - ^34 X 4 = 


X 3 



0) 



Rewriting (9) with b,= elements of fifth column, and a = a , we have 

15 i] lj 



1 


a 





1 









'13 



1 









b„ = = x o 

35 3 , 



V. To obtain.the next matrix a transformation similar to step IV is applied 
to (9): 

1 a. 



'12 





















b 15 - a !3 X 3 

b 25 ' a 23 X 2 =X 2 



Rewriting (10) as was done above with (9), we have 



1 5 12 



VI. Finally, repeating as before, we obtain 
(l 



b 25 =x 2 



b l 5 - a !2 X 2 = X l 



■). 



(9) 



(10) 



(10) 



(11) 



Matrix Inversion: 



Consider (1) with the unit matrix I replacing the column vector b: 



41 






1 














1 





o„ 








1 



45 



48 



d') 



29 



Proceed as in the case of linear equations step la) and b) to obtain: 



/ 


/. 


a i 2 

a 
11 


1 

a a 
ii ii 













a 21 


V ' 


• a 2 4 


1 








\ 


^41 


a 4 2 - • 


* a 44 








1 


and 














/■ 


In 
a ll 


. 


1 
a il 








( 





a 22 - 


-fe) • 


■ -- : ft) 


1 








V 


a 4 2 - 


'"ft) • 


•— w 








1 



(2') 



(3') 



Here again transformations (2*) and (3*) can be used instead of (2) and (3). 
Proceeding as in steps n and III a matrix corresponding to (8) is obtained: 



13 



1 





2 4 



34 



35 



45 



26 



36 



17 




~b„ 



(8') 



48 



Corresponding to (9) we have 




13 






b is _ a i4 


K 


" a i4 


b 4 S 


_a !4 


b 47 


" a !4 


K 





' 5 25- a 24 


b 45 


b _a 

26 24 


_b 4 6 


" a 24 


b 

47 


-a 

24 


b 

48 





b -a A 

35 34 


'K 


b -a A 

36 34 


_b 46 


b -a 

37 34 


b 

47 


-a 

34 


b 

48 



(9") 



where again we can write b for the elements of columns five to eight and a for a . 
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Finally corresponding to (11), the following is obtained: 



1 b-a, o b b-a, b„. b,-a, b„ 

I 15 12 25 16 12 26 17 12 27 

Rewriting the entire matrix as before, the inverse results: 

- • b. 



38 



48 



111 








0l4 


b is 





1 








^25 








1 





b 3s 


°« 








!«. 


b 45 



b 18- a i2 b 28 (11') 



(ID 



Remarks: 

a) Check sums: An entire matrix is augmented by an additional column whose 

elements are the sum of the elements of the respective row: i.e., S = 2a,., where 
S represents the check sum of the i th row. After transformation of the matrix, the 
sum of the transformed row elements is equal to the transformed element of the 
sum; i.e. , 2 T (a 4 ) = TSj , due to linearity. 

/3) Linear equations with the same matrix of coefficients but different constant 

terms (b 's) can be treated simultaneously by forming matrix (1) augmented by several 
column vectors of constant terms. Alternately we can solve for the inverse matrix A -1 
and then perform matrix -vector multiplications x = A -1 b for different sets of column 
vectors b. 

y) The "back-solution" , i.e., the reduction of the triangular matrix (8) to the 

unit matrix, can be eliminated by the use of the following composite matrix: 



-I 











If we apply transformations (2) and (3) [here we cannot use (2*) and (3*)], we obtain 
the values of the unknowns x in the column of the original O-vector located under the 
b-vector, and the inverse matrix A 1 will appear in the original 0-matrix located 
under the I-matrix. (Reference: Verzuh, F. M. , The Solution of Simultaneous Linear 
Equations with the Aid of the 602 Calculating Punch, Mathematical Tables and Other 
Aids to Computation, Vol. 3, No. 27, July 1949). 

The back solution can be avoided without the necessity of having to add the 
-I and matrices by means of the following scheme: Following the first reduction 
the first row of the reduced matrix (3), instead of being discarded, is left in that matrix, 
but is relocated as its last row. (The first row is put below the last row for the pur- 
pose of uniformity in machine operations). The leading column is removed, as before. 
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The reduction is then performed again, this time with the original second row used as 
the first row. After this second reduction is complete, the original second row 
becomes its new last row, etc. Thus each row, in turn, is used as the leading row 
during one reduction, and following that reduction becomes the last row of the reduced 
matrix. The unknowns x will appear in the columns of the original b vector, and the 
inverse A" 1 will appear in the original I-matrix. 



B. Machine Procedure 



1. Consider the matrix 



a il a i2 a i3 a i4 



b 15 K \l 



a 21 a 22 a 23 a 24 b 25 K 27 

a 3! a 32 a 33 a 34 b 35 b 36 37 

a 41 a 4 2 a 43 a 44 b 4 5 b 46 47 



"V 



°ie °X9 °i,io 

^s 29 2)10 

°3 8 1 S9 Sjlo 

48 49 1 4 } 1 



~ ^ 1,11 

-s. 



2,11 



-2 



3,11 



-S 



Elements in columns 1-4 represent the matrix of the set of four simultaneous equations. 
Columns 1-4, 5 and 1-4, 6 are the augmented matrices representing the simultaneous 
equations with two sets of constant values . Both augmented matrices may be solved 
simultaneously. Elements in columns 7-10 represent the unit matrix which will become 
the inverse of the original matrix when the same transformations that will transform the 
original matrix to the unit matrix are applied to elements in columns 7-10. Elements in 
column eleven represent the negative check sums. The matrix above can be augmented 
by the negative unit matrix and the zero matrix as previously described (see previous 
section if) for the purpose of avoiding the back solution. 

As many sets of equations may be solved as desired. Each set is identified by a different 
equation number. Every element of each matrix to be solved is punched on a separate 
card. 

For matrices with less then ten rows, one card column is sufficient for the row number. 
Two columns have been set aside here to indicate the procedure for larger matrices. 

Punch the cards in the following manner (assuming each element to be an eight-digit 
number): 
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Col. 1-2 


Col. 3- 


•4 


Col. 5-6 


Col. 7-14 


Eq.# 


Row # 


i 


Col.#j 


Matrix eleim 


01 


01 




01 


± 
.xxxxxxxx 


01 


01 




02 




01 


01 




11 




01 


02 




01 




01 


02 




02 





2. The reduction of the matrix involves the repetitive applications of transformations 
(2) and (3). This can be done on the machines by repeating the following sequence of 
four main steps: 

a. Sort the cards in sequence by columns (j) (i. e. , sort on row# i first, then on 
column*], finally on equation #). 

b. Using row 1 cards as master cards, gang punch its elements into cards of the 
corresponding column, i. e. , 

Gang punch 



1 



.21 
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c. Sort the cards in sequence by rows (i) (i.e. , sort on column#j first, then on row 
#i, finally on equation#). 

d. Using column 1 cards as master cards, obtain on the calculating punches 602-A 
and 604 

a i a 



ii 



where a n and a n will appear on master cards. 
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Remarks: 

a) Master cards will have to be X-punched preparatory to step b and d above." 

fi) After each reduction of the matrix to one of a lower order, the new elements 

a can be reproduced into a new set of cards into the same field where the original a 's 
were punched. Thus none of the control panels used in steps a to d have to be altered. 
Alternately reading and punching fields can be re -wired on the control panels and new 
elements can be punched successively to the right of the card thus reducing the number 
of cards used. After each reduction step the cards corresponding to the first row and 
column are set aside and the remaining cards represent the matrix of the next lower 
order. Only the first row of the composite matrix (-1, 0) need be used in the first 
reduction step; for the second reduction, the second row is added, etc. 

y) Check sums should be verified after each reduction step when the cards are in 

sequence of rows. Summing the elements in each row on the 405 should give a result close 
to zero since the negative of the check sum was punched in the original cards. 



6) Step d is first performed on the 602 -A and is then verified on the 604 

AB 
by means of the DPBC device. Note that operation of the type ---- can be done on one 

program step on the 602-A. The linear term a ( . can be read into the counter group, where 



a a 
ii ij 



a 
ii 



will be developed at read time; a should be used as the dividend since a appears 
only on master cards and cannot be read directly into the dividend counter. The 
604 operation should be programmed to divide first and then to multiply so that the 
calculation is done in exactly the same manner as on the 602-A, thus avoiding any 
possible discrepancy in the unit position of the result. 



An alternate method, avoiding X-punching, is suggested by William D. Bell in his article 
on " Punched Card Techniques for the Solution of Simultaneous Equations and other Matrix 
Operations, " Proceedings, Scientific Computation Forum, pp. 28-31, New York: IBM, 
1948. 
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MATRIX AND VECTOR ALGEBRA 

Eric V. Hankam 
Watson Scientific Computing Laboratory 



A. Mathematical Preliminaries 

Definition: A matrix is a rectangular array of elements: 



A = (a i4 ) = 



a ,i a, 

11 12 
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ml 



In particular, if n = m, we refer to it as a square matrix. 

An n-tuple of elements ( a, a, . . .a ) is called a vector: thus a 1 x n 
matrix [ a 1 a g ...a ] is called a row vector, and an nxl matrix f a i~l is 



called a column vector. 

Matrix Algebra 

1. Addition of Matrices: 
Consider two matrices, A = ( a r ) and B = ( b ). Their sum A + B is defined 
as that matrix C ( c ) obtained by adding corresponding elements: 



ml 



b„. • • b, 

11 In 



b • • • b 

ml mn 



a + b 
ii T ii 



a ,+ b , 

ml ml 



In In 



a + b 

mn m n 



i.e., (a u ) + (b 1} ) = ( a|] + b tJ ) = (c u ). 
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2. Scalar Multiplication: 

Consider a matrix A = ( a ) and a constant c. Their scalar product cA is 
defined as that matrix obtained by multiplying every element a tj of A by c: 

c- ( a lf ) = (ca,. ). 

3. Multiplication of Matrices: 

Consider two matrices A = ( a tj ) and B = ( b^). Their product AB is defined 
as that matrix C = ( c ) obtained in the following manner; 

n 
k = l 

We see from the definition that the product AB is defined only if the number of 
columns of A equals the number of rows in B. Furthermore, note that 
AB / BA. The multiplication of two vectors is frequently called their "inner 
product". In accordance with the above definition, the inner product is the 
sum of products of the elements of the two vectors . Thus to obtain the product 
of two matrices A and B we have to find the inner products of every row 
vector of A with every column vector of B. The product of a matrix and a 
column vector arises frequently in the theory of linear equations. 

4. Inversion of Matrices: 

The inverse A~ of a matrix A is defined as that matrix which satisfies 
the following relationship: 

AA" 1 = A _1 A =1, 



where I is the unit matrix 1 

9 l 

6 • • -1 

An inverse is defined only for non-singular matrices, i.e. matrices with 
non-zero determinants. 
B. Machine Procedure 

1. Addition of Matrices: 

This can be done on the calculating punches or the accounting machine, 
a. Each matrix element a 1 is punched on a separate card with identification 
i and j. Similarly each b is punched on a separate card. Cards are then 
arranged in sequence of i, j in such a way that corresponding a and b^ 
cards follow each other (sorter, collator). The operation a^ + b tJ = c 4j 
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can then be performed on any of the calculating punches. The result c can 
be punched on the b cards, or, alternately, trailer cards can be merged 
behind each pair of a and b cards, and c can be punched on the corre- 
sponding trailer card. If the operation is performed on an accounting 
machine the sum a + b can be printed or summary punched on both. The 
jnethod can readily be extended to the addition of several matrices, 

A + B+...+ H= (a +b„ + .. . +h ) 
At times it may be convenient to punch one pair of elements a and the 
corresponding b i . on the same card prior to their addition, 
b. Several elements a., are punched on one card. If the order of the matrix 
is not too large, an entire row (or column) of the matrix can be punched on 
one card. In general, the number of elements to be punched on a single 
card will depend on the capacity of the machine to be used for the addition as 
well as the order of the matrix and the size of its elements. For maximum 
efficiency, the number of elements punched on each card should be such that 
they can all be summed in one run without any change in the control panel 
wiring. 

2. Scalar Multiplication: 

This is an extremely simple operation that can be performed on any calculating 
punch. The scalar c can be stored from a master card, and used as a group 
factor on every matrix element a . The only consideration is to determine 
how many such multiplications can be performed per card. This again will 
depend on the machine, the order of the matrix and the size of the a . 's. If 
only one matrix is involved, or if the scalar c is the same for all matrices, 
c can be digit emitted and does not have to be read into the machine on a master 
card. Incidentally, it is implicitly assumed throughout this discussion that the 
various calculations are performed simultaneously on several sets of different 
matrices. This, of course, will require a further code identifying each card as 
belonging to a given set. 

3. Multiplication of Matrices: 

I. a. Each element a is punched on a separate card with identification i and j. 
Similarly each b r is punched on a separate card. We then proceed to re- 
produce the a 4 .'s (in sequence of rows) into the b i 's (in sequence of columns). 
Since each row of the A matrix has to be multiplied by each column of the B 
matrix we need as many duplicates of the B matrix as we have rows in 
the A matrix. Thus any given row of the A matrix will successively 
be paired off, i. e. reproduced on every column of the B matrix. Each 
card will contain a pair of elements a and b , . Trailer cards (one 
for each c ) are then inserted behind each group of cards containing 
elements of one row of A and one column of B. The product elements 
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c = S a .. b are obtained on any calculating punch, and are punched 
into the trailer cards, 
b. Several elements b are punched on one card. If the order of the matrix 
is not too large, an entire row of the B matrix can be punched on one 
card. Then each row of the A matrix has to be reproduced only once 
into a given set of B cards. This requires as many duplicates of the 
B matrix as we have rows in A so that any given row of A will appear 
paired off with a set of B cards containing all the corresponding column 
elements. Then, upon merging of trailer cards, several Sa ik b 
(for different b's) can be obtained simultaneously on the calculating punch. 
It is desirable to select the number of b 's punched per card in such a 
way that all the sums of products can be obtained in one run without any 
change in the control panel wiring. A combination of methods La. 
and I. b. is often indicated. This is, in general, the most efficient 
manner for matrix multiplication. On the other hand punching as 
many b 's on one card as the card will hold, will reduce the number 
of cards, but will require a rewiring of reading and punching fields of 
the control panel. 

If the order of the matrices is sufficiently small we can not only 
punch an entire row of the B matrix per card but also reproduce all 
rows of the A matrix into corresponding cards of the B matrix (in 
sequence of columns). Then all c 's can be obtained on a single trailer 
card. In the extreme case (e.g. .matrices of order 2) all elements of 
both matrices can be punched on one card. 

Method I. b. essentially suggests to pair off one row of the A 
matrix with several columns of the B matrix. Conversely, we can 
reproduce several rows of the A matrix into cards containing corre- 
sponding elements of a single column of the B matrix. Which one of 
the two procedures is adopted depends on the rectangular shape of the 
matrices. For example, for an A matrix with more columns than 
rows, i.e. , long rows and short columns, and a B matrix with more 
rows than columns, it will be preferable to pair off one row of A with 
several columns of B. 

Finally, it should be remarked that it is possible to obtain sums of 
products on an accounting machine by means of "progressive digiting". 
This method can be useful if the matrix elements are of small size. 
II. a. and b. This is essentially the same method as I. a. and I.b. Instead 

of reproducing elements a . of A cards into B cards, A cards can be 

interspersed with B cards on the collator. The factors a and b are 

ij ij 

then read from alternate cards into the calculating punch. 

In order to avoid duplicating B cards, (i.e. , use only one set of A 
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cards and B cards), it may at times be convenient to do several runs on 
the calculating punch, and to follow each operation on the calculating punch 
by a collator run; each time the A cards and B cards are interspersed in 
a different manner until each row of A has been paired off with each column 
of B. 
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PART 2 



IBM Type 604 Electronic Calculating Punch 
and 
Type 602 -A Calculating Punch 
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NOTES ON THE IBM TYPE 604 ELECTRONIC CALCULATING PUNCH 
AND 602-A CALCULATING PUNCH 

Eric V. Hankam 
Watson Scientific Computing Laboratory 

SYMBOLIC NOTATION 

Pj = i th program step 

P = read cycle 



R 


= 


read ; read cycle 


P 


= 


punch ; punch cycle 


RI 


= 


read into 


RO 


= 


read out 


RE 


= 


reset 


RR 


= 


read out and reset 


— 


= 


read into 


=>• 


= 


yields (implies) 


€ 


= 


in (belongs to) 


MP 


= 


multiplier 


MC 


= 


multiplicand 


MPL 


= 


multiply 


DIV 


= 


divide 


N 


= 


normal 


T 


= 


transferred 


X-p.u. 


= 


X-pickup 


D-p.u. 


= 


digit pickup 


I-p.u. 


= 


immediate pickup 


Pil-Sel 


= 


pilot selector 
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604 

F = factor storage 

G = general storage 

MQ = multiplier-quotient storage 

EC = electronic counter 

Ri = read units into 

Ro = read units out of 

Em = emitter control 

B.T. = balance test for step suppression 

Sup = suppress 

Sup + = suppress on plus balance 

Sup - = suppress on minus balance 

C. C. = card cycle 

Pil-Sel = pilot selector 

Pch-Sel = punch selector 

Cal-Sel = calculator selector 

602 -A 

[ ] = counter 

( ) = storage unit 

Co-Sel = co-selector 

B. T. = balance test 
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a) 604 



Multiplication: 

In multiplication, the MC cannot be stored in either EC or MQ; i.e., a 
factor that has been developed in EC must be transferred to another storage unit 
before multiplication by it can be started. It is not possible to square a number 
by reading it out (as MC) from MQ. Similar statements hold for division. It 
is possible, however, to read a single -digit MC directly from the emitter control 
and impulse MPL at the same step. It is not possible to shift during multiplication. 
(Multiplication} is performed internally by over and over ( addition 1 
| Division j ~ (subtraction). 

Pilot Selectors and Coupling Exits: 

Thel-p.u. and coupling exit of a P-Sel are common hubs. (They can be 
regarded as either outlets or inlets) . 

Unfinshed Program: 

An "unfinished program" should be wired to "Stop 1 st Reading" at all times 
as a precautionary measure. If there has been an unfinished program condition, 
punching is suppressed from the calculating units, but the gang punching is not 
affected. 

Storage Assignments: 

It is not possible to select storage assignments during calculate time on a 
program step. 

By an 8-6 assignment any 3-position storage unit can be associated with any 
5-position storage unit. An 8-6 or 6-4 assignment disassociates the storage 
unit involved from its sign control unit. If used independently (i.e. .without an 
associated storage unit) to RI or RO a number the sign will not travel with it. 

Storage Units: 

All storage units except EC read true figures with proper sign at all times 
(internally). In EC a true figure is represented by a complement and conversely. 
Avoiding overlap, two distinct numbers can be stored in EC provided that their 
signs are alike. It is possible to read out from one storage unit into several 
other storage units on the same program step. If a positive number is transferred 
from one storage unit into another storage unit and EC on the same step, it will 
enter the storage unit as a complement. This can be used as a method for con- 
version of numbers to complements. A negative number, however, will read into 
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the storage unit and EC as a negative number. Two distinct numbers with any 
possible signs can be read into a single storage unit (as positive numbers) at read 
time (avoiding overlap). Their signs can be "stored" by means of Cal-Sel. 

It is possible to read from any general storage unit, MQ, or EC into 
F at [puncnl time b y wirin S all positions of, say, G2-F2 on the double-panel 
board. G2 RO and F2 RI are wired from C.C. This may be particularly 
useful on machines without factor storage read-in during calculate time. 

Impulsing Storage Units: 

To RI from an X-card only, a storage unit can be impulsed directly by the 
X from the reading brushes. In this case, column splits must be used for digits 
over-punched with an X. An alternate method that does not require the use of 
column splits, is to energize the I-p.u. of a Pil-Sel with the X, and to take a 
C.C. pulse through the T-side of the selector to the storage RI hubs. A C. C. 
pulse should never be taken through the N-side of a selector (energized by the 
I-p.u.) because the early part of the C.C. pulse would get through the selector 
every time before it is picked up (unless a Y-punch is used to energize the Pil- 
Sel). 

Storage RO hubs (2-panel board) accept impulse for the same length of 
time as C. C. (12. 7 to 9. 7). They must receive an impulse starting before X- 
time and extending until past 9 -time. 

It is sufficient to impulse any storage RI with an X-pulse (11 to 11. 3 time) 
even though the storage RI hubs are inlets from 14. 9 to 11. 7 time. But with the 
exception of the units position, the storage unit will not register any digit over- 
punched with an X unless the X is filtered out by a column split. 

Emitter Control: 

If more than one digit is to be entered into a storage unit by means of the 
emitter control, the digits must be entered first into EC, since only one digit 
at a time (on a single program step) can be entered, and a storage unit would 
reset the i th digit upon entering the (i +l) th 

Double Punch Blank Column: 

To suppress DPBC for cards with an X in column i wire "DPBC" to "Stop 
2 nd Rdg" through a Pil-Sel controlled by X-p.u. from 2 nd reading brushes. The 
delayed pick-up is necessary because the DPBC hubs emit an impulse on the cycle 
following the DPBC operation at the second reading station. 



44 



Checking: 



To check that two card columns are both either X-punched or NX-punched, 



wire: 



PIL-SEL 
X-P.U. 



COLUMN 
SPLITS 




STOP 2ND READING 
O O •* 



C.C. 



T *xf T 

N «T T» N 1 

I 2 

PIL-SEL 



COL i COLj 
2ND READING 



By wiring to "Stop 2 nd Reading" the error card will be the last one in 
the stacker when the machine stops. No error light appears. 



Shifting: 



It is possible to shift to the "left" directly from one storage unit to another 
without going through EC; e.g., 



G2 



G4 




G2 RO; G4 Rl; Ri 4TH 



However, it is not feasible to shift into the other direction without going through 
EC. 
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COMMON ENTRY CHANNEL (R. I.) 

T 



-T-r 



ENTRY 



GS2 



EXIT 



-T-t 



I I I 
I I I 



54321 54321 



ENTRY 



MQ 



EXIT 



54321 54321 



COMMON EXIT CHANNEL (R.O.I 



COUNTER 
ENTRY 



13 12 II 10 9 8 7 6 5 4 3 2 I 

COUNTER 

13 12 II 10 9 8 7 6 5 4 3 2 I 



COUNTER 
. EXIT 



13 12 II 10 9 



wu 



I 



"ttttttttrtttt 

13 12 II 10 9 S 7 6 S 4 3 2 I 



8 7 6 5 432 I 



Referring to the above figure, in the normal position of the shift unit 
the first five exit positions of the counter are lost in a transfer out of the 
counter (i.e., positions 1, 2, . . . , 5. ) This position of the shift actually 
corresponds to Read Units out of 6 th . To effect a normal units to units 
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transfer out of the counter it is necessary to shift the units to the farthest 
position to the left. This will automatically occur when a counter read-out 
is programmed. 



Therefore: 




Normal R. 0. of counter = 


Ro 6 th 


Ri 2 nd = 


Ro5 th 


Ri 3 rd = 


Ro4 th 


Ri 4 th = 


Ro 3 rd 


Ri 5 tb = 


Ro 2 " i 


Ri 6 th = 


Ro Normal 



Note: Corresponding hubs are common, i.e. , Ri 2 td and Ro 5 th are common hubs, 
etc. 
Timing: 

Calculate time may last beyond "Y" time; therefore it is seldom safe to 
use a Y-punch. Since "Y" time precedes card cycle time, it can be used to select 
the C.C. impulse. 
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10 X 10 MULTIPLICATION 
Examples: 604 



A 
5 5 



B 
5 5 



R : a, - 



a„ — 



b, - 



b„ - 



F2 
MQ 
F4 
Fl,3 



p 9 '• 



p 10 : 



a 2 b 2 : Fl,3 RO 



: a 2 b 2 -Gl,3: EC RR 
: a 2 b- EC : Gl,3 RO 



a 2 b i : F4 RO 



a,- MQ : F2 RO 
aj b 2 : Fl,3 RO 

round : 1/2 Adj 

a' b'-Gl, 3 : EC RR 
a' b'-EC : Gl, 3 RO 



a x b l : F4 RO 



5 
I 1- 



J i a b 



a 2 b, 



a A 



±5 (ROUND) 



MPL + =». a 2 b 2 e EC 

Ro 6th ; Gl, 3 RI (highest 6 positions only) 

ECRI + 

MPL + ^ a 2 b 2 + a 2 b x e EC 

MQ RI 

MPL + =* ^ b 2 + a 2 \ + ^ b 2 = a' b 'e EC 



Ri 5th x xxxxx I xxxxx 
carry ±5 

Ro 6th ; Gl, 3 RI (highest 6 -positions only) 
EC RI + 

MPL + <*. a 2 b 2 + a 2 b x + a 1 b 2 + a ; ^ = AB e EC 



P : AB (10 digits) from EC 
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15 X 15 MULTIPLICATION 
Examples: 604 



R: 



a 

FI.3 



1 

b c 

i i 

F2 MQ 



-+- 



ABC 

i I I 

Gl,3 F4 G4 



P 2 : 



V 



P, :' 



V, 



P s : 



p 8 = 



cA 
C -MQ 

aC 

B - MQ 
Bb 

f round 

A-G4 
M r Gl, 3 



V 



P l2 : 



13" 



P !4 : 



P i 5 ; 



*V 



p 1T : 



Mj- EC 



Ba 
A-MQ 

bA 



M-G2 

2 

M-G1.3 
M- EC 



aA 

10 
xxxxxxxxxx 

EC 



Gl, 3 RO 

G4 RO 

Fl,3 RO 

F4 RO 

F2 RO 

|Adj 

|Adj 

Gl,3 RO 

EC RR 

Gl, 3 RO. 

Fl,3 RO 

G4 RO 

F2 RO 

EC RO 

EC RR 

Gl,3 RO 

Fl,3 RO 

5 
xxxxx 

G2 



Gl, 3 RI 



MPL + 
MQRI 
MPL + 
MQ RI 
MPL + 
Ri 5th 
Ri 5th 
G4 RI 
Ro 6th 

EC RI + 

MPL + 

MQ RI 

MPL + 

G2 RI 

Ro 6th ; Gl, 3 RI 

EC RI + 

MPL + =► 

aA 



►cA e EC 



-cA + aC e EC 



-cA + aC+ Bb = M 1 e EC 



x xxxxxlxxxxx 
carry ±i.o 



(highest 6 digits only) 



•M + Ba e EC 



»M + Ba + bA = Me EC xxxxxxxxxxx 

1 2 



(last 5 positions) 
(highest 6 positions) 

highest 10 positions of product e EC 



h 



5 5 

-I 1 



aB+bA 
cA + aC+Bb 



Note: 



Adjustment of 2. 5 for each of the two 
omitted products cB and bC is necessary. 



H 



cB+bC J_ 1 .] 

cC ! t | 1- 1 

ROUND 5-1-2.5+2.5 = 10 
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16 ± 16 
Examples: 604 



16-9 



'16-9 



R : Gl,2 



p i : a i6- 9 - EC 



p 2 : b 16 -9 - EC 



P : a + b - Gl, 2 

3 I6-9 ' 



P 4 : a+ b 17 -G3 



P 5 : preserve sign of G3 
P 6 ' balance test 

Sup+P, : +1 - G4 
7 (if<o) 

P : -1 - EC 

8 

P 9 : balance test 
P 10 : clear EC 
.Sup+P^: -1 - EC 

Vftr EC 

P i 3 : balance test 
P 14 : clear EC 

P 15 : a + b 17 -EC 

p 16 : clear G3 

p . cancel extraneous 
" '. MQ value 

p ia : a + Vr EC 

Sup"P : "boost" signed 
quantity 



I 
FI.2 

Gl,2 RO ; ECRI + 

G3, MQ RO ; EC RI ± 

EC RO 

EC RO 

EC RO 

B.T. 

Em "1" 



G3.MQ F3,4 



a ± b 17 _ g € EC 



Gl,2 RI 

Ro 4th ; G3 RI 

Ro 4th ; MQ RI 

04 m (to clear G4for suppressed case) 

G4RI =► {°} e G4 

EC RI- : 



a + b 17 ^-l e EC 



Em "1" 

B.T. 
EC RR 
Em "1" 

G4 RO 
B.T. 
EC RR 

G3, MQ RO ; EC RI + ; Ri 4th 

G3 RI 

G3, MQ RO ; EC RI- ; Ri 4 th 

Gl,2 RO ; ECRI+ =*-a + b e EC 

17-9 



; EC RI- => ^ o) e EC 

; EC RI +=*[_?] e EC [-1 £ E only if a + b W _ B = 0] 

[ note: if < 0, a + b,, = ] 



x 00000000 e EC 



i Adj 



x xxxxxxxx 

±5 
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Sup-P, 



31 



34 



a + b 17 _ 13 - MQ 



a + b „ , - EC 

17-13 



a + b - Gl, 2 

12-9 



a + b iM -EC 



a+b i 2 ^ G1 ' 2 



a + b,,^- EC 

12-9 



Vi - EC 



b a-i - EC 
a + b c ,- G4 

5—1 



a + b is-, "* G1 ' 2 



a+b X3^ - EC 



cancel ± 5 



ja+b^ -EC. 



a + b„ - Gl, 2 

17-13 ' 

a + b — EC 

«■ T "17-13 "^ 



EC RO 
G3, MQ RO 
EC RR 
Gl, 2 RO 
EC RR 
Gl,2 RO 

Fl,2 RO 
F3,4 RO 
EC RO 
EC RR 
Gl, 2 RO 
|Adj 
Gl, 2 RO 
Gl, 2 RO 
G3,MQRO 
Gl, 2 RO 



Ro 5th 
ECRI- 
G1,2RI 

EC RI + 
G1,2RI 
ECRI + 

ECRI + 
ECRI ± =» 
G4RI (for 
Ro 6th 
ECRI- 
Ri 4th 
ECRI + 
EC RI + 
Gl, 2 RI 
ECRI + 



; MQ RI 

; Ri 5th subtract out 00000 xxxx 

a+b 



12-9 



; Ri 4th 



; Ri 6th xxxx 00000000 e EC 



• a + b x xxxx xxxxxxxx e EC 
"-carry 

punching a + b from G4 with proper 

sign) 
; Gl, 2 RI 

xxxxx'xxx 

±5 (from P ) 

*5 



• a + b l3 _g e EC xxxxxxxx 



a ± b,„ . from EC 



; a± b from G4 



Ri 4th 

Ri 5th=*-a + b M 
(use sign of G4) 



,1,2 
e EC xxxxxxxxxxxx 



Remarks: 

1) P 5 serves to " store " the sign of a + b 17 (in G3) in MQ for the transfer 
of a + b at P . The extraneous value in MQ is subtracted out at P . 

G3 is cleared at P 16 and, although "coupled" to the MQ RO, does not 
interfere with any subsequent MQ RO. 

2) P is required to increase the number in EC sufficiently so that a + b 
will have the proper sign when it is read into G4 at P 28 . P al subtracts 
out the ±5 added at P,„ . 
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(continued on next page) 
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3) a 16-9 + b 16 ^ = is treated as a special case. P lg and 1^ are suppressed 
whenever a, + b, = a + b =0. 

16-9 16-9 16-9 



4) TO test whether a number = in EC: 

B. T. and if < 0, emit 1 into a storage unit, say G4. 



EC B.T. G4 

<o - +i 



Subtract -1 to the number in EC . 

B.T. ; clear EC ; and if <0, emit -1 into EC. 



Read G4 -EC: 



EC BT EC 




>° } + 
J 




<o - -l 




G4 EC 




>0 case : + = 





case : - 1 = 


-1 


<0 case : +1 - 1 = 






B. T. again : Only zero case will give negative balance. (The zero case is 
characterized by + balance after first and - balance after second B. T. ). 
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b ) 602 -A 

Multiplication: 

It is not possible to perform two multiplications (or two divisions, or one 
multiplication and one division) on successive program steps. At least one pro- 
gram step must intervene between the multiplications. 

When squaring a number it is possible to put the number in (1R) and use it both 
as a MP and MC; 

It is possible to wire a constant MC directly from a digit emitter on the 
step at which multiplication is performed. 

Depending on the digits of the MP, several cycles are required during a 
multiplication (similarly division) . The first multiplication cycle occurs during 
the program step on which multiplication is impulsed. That program step can be 
regarded as "repeated" for each of the remaining multiplication cycles. 

A Co-Sel (or Pil-Sel; I-p.u.) picked up by a program couple or program exit 
(of the same step on which MPL is impulsed) will remain transferred throughout 
the entire multiplication. 

If the MC is in a counter, do not reset the counter during the multiplication 
step. It can be reset on the step following the multiplication. 

A counter should not be reset on a program step during which multiplication 
takes place, even though the counter is not involved in the multiplication. 

Couple adjacent counters for the development of a product or quotient. 

The digit emitter is active on each multiplication cycle; it is inactive for a 
zero multiplier. 

Therefore a Pil-Sel energized by a digit on a program step on which a multi- 
plication is performed will pick up on the second multiplication cycle, but will 
not pick up at all if the multiplier is zero. A Pil-Sel needed for the extension of 
program steps should not be picked up by digit impulses at P if a multiplication 
takes place at P . 

12 

Balance Test: 

Once a number is punched converted as a result of a B. T. , all following 
results for the same card will also be converted until a new B. T. takes place. 

A B. T. unit can be used more than once by split wiring. 

More than one B. T. can be made by two distinct numbers being punched 
from the same storage unit on the same program step. 

Zero (in a counter) is regarded as a positive quantity during a balance test 
operation on the 602, 602-A and 604. 

Punch Interlock: 

Wiring P i to "punch interlock" will hold up P £ until the previous card has been 
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completely punched. 

Punch Interlock cannot be used to delay a program step until punching 
(impulsed on a previous program step) on the same card has been completed. It 
also cannot delay the read cycle of the following card. 

Punch Emitter: 

Do not split wire two hubs from the (external) punch emitter into a single 
punching position. 

The punch emitter is not timed; it emits as soon as the card moves to the 
particular punching position. 
Counter: 

A true figure appears in a counter (internally) as a 9's complement. A 
counter filled with 9's is treated as zero, (holds for 602 and 602-A). 

All counters can be coupled together into a 40-position counter group in 
which 40-digit numbers can be accumulated. 

It is possible to accumulate the product of two numbers in two different 
counter groups at the same time. In the case mentioned in paragraph 2 under 
"Multiplication", a 2 was computed, and +a 2 was accumulated in one counter while 
-a 2 was developed in another counter. 

Card Feed Operation: 

There are two cycles to a card feed. A card passes the reading brushes 
on the first one, the read cycle, and the card behind it passes the control brushes 
on the second one, which is normally program 1. Thus.it is not possible for 
cards to pass the reading brushes on successive cycles. If all program steps 
are skipped for a card ("Read Cycle" wired to "Read" hub, usually through a 
selector) the read cycle hubs will become active again on the cycle after that 
card passes the reading brushes. However, the following card does not pass the 
reading brushes on this cycle; it is passing the control brushes on the second 
half of the previous card feed operation. 

Transfer of Numbers: 

A number must not be transferred from (6), (7) or (8) on the same step 
when "punching" is impulsed from these units. Transferring a number while 
it is punched should also be avoided. It may cause the occasional skipping of a 
column to be punched. 

It is not possible to transfer numbers during multiplication (or division) . 

Co-Selectors: 

To pick up the same Co-Sel at several (different) program steps, 
impulse the Co-Sel pickup hubs directly from the program exits. 
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Program Steps: 

Method to extend program steps beyond 12: 



COUPLE EXIT 




CO-SEL P.U. 



o 



' 2 
O 



EXITS 



•*,C (DIGIT EM.) 
ANY DIGIT 



: uu 



12 




CO -SELECTORS 
Pil-Sel.l and Co-Sel 1 will be picked up from P x until Read time. 



2 " 

3 " 

4 " 



2 " " 

3 " " 

4 (( t( 



P " 

25 

P " 

37 

P " 

49 



Pgg will get through Pil-Sel 4 to Read. 

Usually several Co-Sel will be coupled to form a sufficiently large group 
to allow taking through all necessary program exits. 

Now suppose we want to start with program steps 9, 10, 11 (following a 
read cycle); then repeat steps 12, 1, 2, 3, 4; 12, 1, 2, 3, 4; etc. , until a negative 
balance condition arises (during one of these steps, say step 3); and finally end 
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up with steps 5, 6 before going back to Read. A wiring scheme for this is as 
follows: 



pil-sel x p.u. 



READ i 

DROP OUT 



OT 



O N CO-SEL 
C 



CO-SE 



SKIP O 
COUPLE O 

O 

O 
EXITS 

o 



L P.U.^ 



O 
O 

o 
o 
o 



o 
o 

O 
O 



O O 



O 

o 

©- 

o 
o 



-*• c 







RE 


AD 




J , 


fll _READ CYCLES 








O O 








o o 


o 


o 




o o 


o 


o 




o o 


o 


o 




o o 


o 


o 




p ° p !\ 


o 

P 9 


o 

p , 2 





NEGATIVE BALANCE 
'(OF A COUNTER) 



To skip program steps in the "middle" (i. e., between 1 and 12), wire: 



THROUGH SELECTOR 



OR 




P 5 , P 6 , P 7 will be skipped . 

P and P will occur at the same time. 

4 8 



^i€^ p s*« iP 




P. > P= . P, will be skipped. 

5 6 7 

P, and P will not occur at the same time. 

4 8 



Read Drop Out Impulses: 

The R.D. O. (Read Drop Out) impulses are common hubs, and may cause 
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back circuits if two Pil-Sel are dropped out in a different manner. For example, 
if both Pil-Sel are dropped out with a R.D.O. impulse, and if one of the two is 
also dropped out at the end of a program step, the R.D. O. impulse to the other 
Pil-Sel should be wired through the transferred side of a Co-Sel picked up at 
Read time. 

Shifting: 

Impulsing MPL and TENTH hubs on the same step will cause the product 
to be shifted one position to the right. 

Skipping: 

The SKIP OUT hubs accept impulses only during a read cycle. To skip a 
card on a condition arising during calculation, wire the first position to be punched 
from the storage exits, and the SKIP pulse through either the N- or T-side of a 
Pil-Sel (which reflects the condition) whose common point is wired to the punch 
magnets. If the SKIP pulse is wired through the N-side of the Pil-Sel, it should be 
taken through a Co-Sel (picked up after the condition arises) to prevent it from 
getting through before the Pil-Sel is picked up. 

The Pil-Sel will usually be dropped out by a "punch drop out" impulse, and 
"punch interlock" should be impulsed so that the Pil-Sel cannot be activated be- 
fore punching of the preceding card has been completed. 

Also it should be mentioned that the "Pj-Skip" hub is dead. 

Division: 

It is possible to perform an operation of the type ( _ )C on the same step. 

In addition to the standard division wiring, C is wired into a counter which will 
develop (JL\c . MPL must not be impulsed. 

A zero divisor will cause the machine to punch zeros for the quotient. 

If the difference of significant digits of the dividend and divisor exceeds 7, 
the machine will punch an incorrect answer after very lengthy calculation. 

The highest position of [l] must remain free in division; [l] must always 
be used as the dividend counter. Adjacent counters can be coupled for placement 
of a larger dividend; usually [1,2, 3]. At least one significant digit must appear 
in (1R), the divisor unit. It is possible to use a divisor with more than 8 
digits by placing the low order digits into any storage unit; the dividend can also 
be extended beyond 15 digits by coupling counters adjacent to [l, 2, 3]. 

In any division problem there exists a one-to-one correspondence between 
the units position of (1R) and the units position of [3] in the sense that any signi- 
ficant digit standing in the units position of (1R) must be wired into the units po- 
sition of [3] during division. Digits located n places to the ( l e ^J of the units 

(right) 
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position of (1R) must be wired n places to the j^J of the units position of [3]. 
For example, consider a 14-digit divisor; it must be read into [l, 2, 3,4, 5] in the 
following manner: 



(IR ) (2L) OR ANY OTHER STORAGE UNIT 



AAAAA**\« « " " " " " 

tninn unit 



O 00 00 00 X X X XX X X X xxxxxxoo 

\nITS POSITION OF(IR)** UNITS POSITION OF Q3] 

Thus it can be seen that for an 8-digit divisor located entirely in (1R), 
as in standard division, [l, 2, 3] suffice as the dividend unit; and for a 4-digit 
divisor entered into the highest 4 positions of (1R), only [l, 2] need to be coupled 
as the dividend unit. 

The number of quotient digits (« 8) developed depends on the placement of 
the dividend in the dividend counter group, as well as on the original placement 
of the divisor into (1R). 

Rounding the Quotient in Division:" 

In order to avoid the loss of the quotient position which is used to 
receive the "reset to 5" and then dropped out, it is possible to round the 
quotient in the following way and not sacrifice any positions. (This is, in 
fact, the method of the 602). On the program just preceding the "divide" 
program, wire the four program impulses to 1) storage 1 read-out, 2) x 5, 
3) tenths, 4) dividend counter read -in plus. This will cause half the divisor 
to go into the dividend counter before dividing starts, and will thereby produce 
a rounded quotient directly. 

Division with a Negative Dividend: 1 

If one uses the same method of rounding recommended in the above 
paragraph, and if one has satisfied the rules about not having a divisor so 
small that seven shifts will be unable to test the left end figure of the dividend, 
then the quotient counter receives the correct negative result and this is 
converted to a true figure during the punching. Notice that the machine is 
building the quotient by working the dividend counter up to zero from its original 
negative number and not down to zero from a positive dividend. Nevertheless, 



This process was suggested by Dr. Paul Herget, University of Cincinnati 
Observatory, Cincinnati, Ohio. 
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the machine always leaves a positive remainder. For these reasons, one 
should add 0. 5 x divisor just before the division takes place. 



Punching: 



To punch an X (during a B. T. operation) over the highest position of the 
result, wire: 




PUNCHING 



STORAGE UNIT 



If two numbers are to be punched successively (on different program 
steps) from the same storage unit, an internal interlock is provided to hold 
the second number until completion of punching of the first. For that reason 
the units position of (6), (7) or (8) has to be wired to be punched. It is possible 
to read the first result into, say, (8) at P £ ; impulse it to be punched at P 1+j 
(where j may be 0); read a second result into (8) at P t +j +k (where k / 0); 
and punch the second result at P [+ . +k + , (where 1 may = 0). Trouble may 
arise, however, in the case when the same result is to be punched twice 
(into two different card fields) and it is read into the punch storage unit only 
once . Since there is no second RI between the two punchings, the internal 
interlock will not be effective, and either a sufficient number of program steps 
must intervene to allow for the completion of the first punching, or the same 
number should be read into the punch storage unit (from its previous storage) 
again for the second punching. 

Punching is accomplished by an internal "punch" emitter sending out an 
impulse which is selectively taken to the punch magnets depending upon the number 
standing in the storage position. 

Punch storage delays cannot affect the occurrence of a read cycle; i. e.,care 
should be taken in reading from the card directly into punch storage. Nothing will 
prevent the storage from clearing and reading in a new value before the punching 
of the previous card has been completed. 

A number may be read from one card and punched into succeeding cards 
(gang punching operations). If the master card is wired to skip out without punch- 
ing, the impulse to the punch hub must not come before the master card leaves 
the punch bed. Because of the storage interlock circuits, if this impulse comes 
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before the master card is actually skipped out, an error stop will result. The 
use of a "Read cycles" impulse will cause this. Moving the punch time back to 
the first program will clear up this difficulty. 

Split wiring any punch storage position, except the units position, with an 
X from the punch emitter will cause all columns to the left and right of the split 
wired position (except the units column) to be overpunched with an X. 



Pilot Selectors: 

Even though a Pil-Sel impulsed by an I-p.u. need not be dropped out, it is 
preferable to do so, because a Pil-Sel can easily remain transferred as a result 
of a previous control panel (machine run). 

Use a "9"-pulse to drop out a Pil-Sel at the end of a program step. Do not 
use a program exit to drop out a selector. (The drop out impulse should never come 
between 355 and 50 ). 

When a Pil-Sel is impulsed by the I-p.u. , its couple exit is not active. 

To pick up a Pil-Sel with the sign (X in column i) of a group factor appearing 
on master cards with X in column j, hold the selector throughout the passage of 
detail cards, and drop out before reading the next master card: 



I 



X IN COL. J FROM 
CONTROL BRUSH 

X-P.U. PIL-SEL I 



I DROP OUT 
(READ DROP OUT 



t 



X-P.U. PIL-SEL 2 




DROP OUT 
READ DROP OUT 



X IN COL. I FROM 
READING BRUSH 



Pil-Sel 1 serves as auxiliary selector. 

Pil-Sel 2, if picked up, will remain so from one master card to the next. 
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24 X 24 MULTIPLICATION 
Examples: 602 -A 



a„ — 



a, - 



a. — 



b - 



(IR) 
(2) 
(3) 
(4,5) 



a, 

8 



Z4 

B 8 



h a 3 b 
a 2 b 
a, b 



a 3 b 




: (4, 5) RO 


a 2 - (IR) 




: (2) RO 


a 3 b-(2,6) 




[1,2,...,8]RR 


a 3 b - [1, 2, . . 


-,8] 


(2, 6) RO 


a 2 b 




(4,5) RO 


a x - (IR) 




(3)RO 


a x b 




(4, 5) RO 


ab - (7, 8) 




[1,2,...,8]RR 



[1,2,...,8]RI+ ; MPL 

(IR) RI 

(2, 6) RI ^highest 24 digits only] 

[l, 2, . . . , 8]RI + (into units position] 
[1,2,..., 8] RI+ ; MPL 

(IR) RI 

[1,2,..., 8] RI+ ; MPL 

(7,8)RI ; (7,8)P ; READ 



Note: b is read from (4, 5) into [1,2,..., 8] into the units position through the N-side 
of a Pil-Sel at P 1 and P 5 , and into the 9 th position through the T-side at P . 
The Sel is picked up by the P 7 - couple . 

(l, 2, . . . , 8] is reset to 5 into the 25 th position counting from left, through a Pil- 
Sel picked up at P 3 . 
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32 x 32 MULTIPLICATION 
Examples: 602 -A 



B - (2,3,4) 
a 4 - (1R) 



..r,. 8 8 8 8 

MP I 1 1 1 1 

a, a 2 a 3 a 4 
*" lowest 8 positions of (2) 



B 

32 



HMC 





/ 6 \ f highest 4 position 
s "* v '/ ^-highest 4 
a 2 - (2L,5L) 












' e ' 




t 


32 


4 8 " 

'round 




a t - (5) lowest 8 




v 

PUNCH 


P x 


:a.B : (2,3,4) RO 


[1,2,.. 


. , 8] RI + 


; MPL 


Data card 








p 2 


: a 3 - (1R) : (6) RO 


(1R) RI 






P 

3 


:a 4 B-(6,7,8) : [1,2, . . . ,8] RR 


(6,7,8) 


RI 




P 4 


:a 4 B-[l,2,...,8] :(6,7,8)RO 


[1,2,... 


,8]RI + 




P B 


:a 3 B : (2,3,4) RO 


[1,2,... 


,8] RI + 


; MPL 


P e 


:a 2 -(lR) : (2L) RO;(5L)RO 


(1R) RI 






p 7 


:(H +a 3 )B-(6,7,8) :[l,2 8]RR 


(6,7,8) 


RI 




P B 


:(a 4+ a 3 )B-[l,2,...,8]:(6,7,8)RO 


[1,2,.. 


,8]RI + 




P 9 


:a 2 B : (2,3,4) RO 


[1,2,.. 


,8]RI + 


; MPL 


P 10 


: a t - (1R) : (5L) RO ;(5R)RO 


(1R) RI 







-. a 3 B 



P 13 :a ! B 



P l : AB^(6,7,8) 
Trailer card 



a,B 



a 4 B e [1,2,..., 8] 40 digits 

(highest 32 positions only) 
(into units position) 
a 4 B + a 3 Be [1,2,..., 8] 40 digits 

(highest 32 positions only) 
(into units position) 
a 4 B + a 3 B + a 2 B € [1, 2, . . . , 8] 40 digits 



(highest 32 positions only) 

(into units position) 

(2,3,4) RO ; [1,2, ...,8]RI+ ; MPL =»■ (a„ + a 3 + a 2 + a t ) B = AB e [1,2, ... ,8]40 digits 

READ (Data cards only) 

[1,2,...,8]RR ;(6,7,8)RI ;(6,7,8)P (highest 32 positions only) 

READ 



P u :(a 4+ a 3+ a 2 )B-(6,7,8):[l,2,...,8]RR ; (6,7,8) RI 

P 12 : (a 4+ a 3 +a 2 )B-[l,2 8]: (6,7,8) RO ; [1,2, . . . ,8]RI + 



Remarks: 
a) 

b) 



Due to lack of card space the product is punched into trailer cards inserted 
alternately behind each data card containing the factors (64 card columns) . 

There can be no carries when combining partial products because the entire 
MC is always used; e. g.,a 4 B + a 3 B = (a 4 + a 3 ) B constitutes a 16 x 32 = 48 
digit multiplication. 

(continued on next page) 
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c) Rounding is accomplished by "resetting to 5" the 8 "position (counting from 
the right) of [l, 2, . . . , 8] through a Co-Sel picked up at P u . 

d) The highest 4 positions of (5L) and the highest 4 positions of (2L) must be 
taken through a Co-Sel picked up at P 6 . The lowest two positions of (5L) 
must be taken through a Co-Sel picked up at P . 

e) (6) RO and (7, 8) RO can be coupled; likewise their RI. 

f) (2L) RO must be kept separate from (2R, 3,4) RO to avoid back circuits; in 
addition when reading out (6, 7, 8) - [l, 2, . . . , 8] at P 4 , P 8 and P 12 , the highest 
4 positions of this 32-digit number (located in (6)) have been taken through a Co- 
Sel picked up at P 4 , P a , and P 12 . (For method to pick up the same Co-Sel 

at more than one program step, see previous section on 602-A Notes, section 
headed "Co-Selectors".) 

g) Since ^ and P (data card) are almost identical, only the impulse back to 
READ has to be wired through a Selector but P 1 for data cards and P y for 
trailer cards are entirely different, and the distinction has to be made through 
a selector. 
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Ji = 22 DIVISION 
22 

Examples: 602-A 



P o :a-[l,2,...,6l 
b-(lR,2, 1L) 



29 
Oxx. . .x 

XXXXXXXX XXXXXXXXXXXX XX 

(1R) (2) (1L) lowest two positions 

;(1R,2,1L) RO ; [l, 2, . . . ,6] RI- ; [7,8] RI + ; DIV «* (h/ 

' '22-15 



first 8 digits of quotient 

e [7,8] 



P 2 :g) -(8) 

v '22-1 



P 

4 

P 



P 



22-15 

Rem, -(3,4) 



: [7,8] RR ; (8) RI 

:[1,2,...,6]RR ;(3,4)RI 

Rem, -[i,2,..., 61: (3,4) RO ;[l,2, . . . ,6] RI+ 

:(lR,2,lL)RO ;[1,2,...,6]RI- 

: [7,8] RR ; (8) RI 



(8)P 



Rem 



[7, 8] RI + 
(8)P 



(punch first 8 digits of quotient) 
(lowest 22 positions of [l, 2, . . . , 6] 
(shift 7 places to the left) 

(punch second 7 digits of quotient) 



Rem 2 -(3,4) :[l,2 6]RR 

Rem 2 -[1,2,..., 6]: (3,4) RO 
52HL2 :(lR,2,lL)RO 

: [7,8] RR 



GL - 



(8) 



clear Rem 



(8)RI 
:[1,2,...,6]RR ;READ 



(3,4) RI 

[1,2,...,6]RI + 
[l,2,...,6lRI- ;[7,8]RI+ ; DIV =*(|) 6 [7,8] 
I ; (8) P 



(lowest 22 positions of [l, 2, . . . , 6] 
(shift 7 places to the left) 



(punch third 7 digits of quotient) 



Note: if rounding is to be included, the case of a partial quotient consisting of all 9's 
must be considered, since in that case a carry should go over into the units 
position of the preceding partial quotient. 
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29 



= 29 DIVISION 



22 
Examples: 602-A 



P o : a - [1,2,..., 6] 
b -(1R,2,1L) 



i • b 
\ :(e) ~< 8 > 

* \ u /29-22 



29 

o£crrTx 

xxxxxxxx 
(1R) 



: (1R,2,1L) RO 
: [7,8] RR 



Re mi -(3,4) :[l,2,...,6]RR 
Ren^ - [1,2 6]:(3,4)RO 

™?™X : (lR,2,lL)RO 

(8) : [7,8] RR 



21-15 



P 7 
P„ 



Rem 2 -(3,4) : [1,2, . . . ,6]RR 

Rem, -[l,2,...,6]:(3,4)RO 



Rem. 



(l)-< 8 > 



(1R,2,1L) RO 
[7,8] RR 



Rem 3 -(3,4) :[l,2, . . . ,6]RR 

Rem 3 -[1,2,...,6]: (3,4) RO 

: (1R,2,1L) RO 



Rem, 



(1L) lowest two positions 

-first 8 digits of quotient 
[1,2,. ; .,6]RI- ;[7,8]RI+ ; DIV =*(£) ) e [7,8] 



(8) RI 



; (8)P 



(punch first 8 digits of quotient) 



(3, 4) RI (lowest 22 positions of [l, 2, . . . , 6] 

[1, 2, . . . , 6] RI+ (shift 7 places to the left) 

[1,2,...,6]RI- ;[7,8]RI+ ; DIV =*(f ) g [7,8] 

S ' 21-15 

(8) RI ; (8) P (punch second 7 digits of quotient) 



(3,4) RI 
[1,2,...,6]RI+ 



(lowest 22 positions of [l, 2, . . . , 6] 
(shift 7 places to the left) 



[1,2,...,6]RI- ;[7,8]RI+ ; DIV =».(|) £ [ 7j 8] 

* ' 14-8 

( 8 ) M ; (8) P (punch third 7 digits of quotient) 

(3j 4 ) RI (lowest 22 positions of [l, 2, . . . , 6] 

[l, 2, . . . , 6] RI + (shift 7 places to the left) 

[1,2,...,6]RI- ;[7,8]RI+ ; DIV =*(|) e [7,8] 



clear Rem^ 



: [7,8] RR ; (8)RI 

:[1,2,...,6]RR ;READ 



; (8) P 



(punch fourth 7 digits of quotient) 



Note: If rounding is to be included, the case of a partial quotient consisting of all 
9's must be considered, since in that case a carry should go over into the units 
position of the preceding partial quotient. 

Remark: In this problem all 80 columns of the card are punched (29+22+29 = 80). 
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2 2 

2 x; S y; 2 xy; S x ; S y 
Examples: 602-A 



Introductory Description 



£ x will be developed in [l] 
Sy " " " " [2] 



S xy 



S x " " 



s y 



[3,4] 
[5,6] 
[7,8] 



All calculations are to be performed on detail cards. 

Each sum is developed progressively in its respective counter group. 

The final sums are cleared and punched on a trailer card. 

The program step required for reading out the final sums must be skipped for all detail 

cards. The program steps required for the calculations can be skipped for the trailer 

card but need not be if the trailer card is blank in the detail card fields containing x and 

y. 

Outline of Program 



p 

o 


: x 


[1] 


and (1R) 








y 


[2] 


and (2R) 






P x 


: xy 




: (2R) RO 


[3,4] RI + 


MPL ==► S xy 6 [3,4] 




2 
X 




: (1R) RO 


[5,6]RI + 


=► S x 2 € [5,6] 


P 2 


: y-OR) 




: (2R) RO 


(1R) RI 




P 3 


2 

:y 




: (2R) RO 
READ (detail c 


[7,8] RI + 
,ards only) 


MPL =** S y 2 e [7,8] 


P 4 


: punch 




: [l]RR 


, [2] RR 


(6) RI ; (6) P 








[3,4] RR 


[5,6] RR 


[7,8] RR 








(7,8) RI 


(7,8) P 


READ 
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8 DIFFERENCES 

(11 digits) 
Examples: 602 -A 



P 3 

p. 



A - 



7 i 

: A - 

i-7 



6 1-1 

: A - 

1-7 



A - 



P 5 : A 



5 1 

P. : A 

6 1-5 



4 1-1 

P 7 : A . 

' 1-5 

4 1 

P„ : A ■ 



p 

9 


: Y' 1 

i-4 


P 10 


1-3 


P u 


1-3 


P !2 


2 1 

: A 

1-2 


P 13 


1-2 


P,. 


:V 



7 1 

: A 



: f. - 



[3,4,5] 
[6,8] 
[1,2] 
(8) 

[3,4,5] 
(7); (8) 
[1,2] 

[3,4,5] 
(6);(8) 

[1,2] 

[3,4,5] 

(5);(8) 

[1,2] 

[3,4,5] 
(4);(8) 

[1,2] 

[3,4,5] 
(3) ; (8) 

[1,2] 

[3,4,5] 
(2); (8) 

[1,2] 

[3,4,5] 
(1) ; (8) 

[1,2] 

[3,4,5] 



[3,4,5] RI 
[6,8] RI + 
[3,4,5] RR 

[1,2] RO 

(7)RI 

(6)RO 

: [1,2] RO 
:(6)RI 

: (5) RO 

:[1,2] RO 
:(5)RI 

: (4)RO 

:[1,2] RO 
:(4)RI 

: (3) RO 

:[l,2]RO 
:(3)RI 

: (2) RO 

:[1,2] RO 
:(2)RI 

: (1) RO 

: [1,2] RR 

(1)RI 

(7) RO 

[6,8] RR 
READ 



i. I 

f i " [ f ,-i + A . 



A .-J = A e t 3 > 4 > 5 ] 



; [1,2] ri + 

; (8) RI ; (8) P 

; [3,4,5] RI- 
; (8) RI ; (8) P 

;[1,2] RI + 

; [3,4,5] RI- 
;.(8)RI ;(8)P 

;[1,2]RI + 

; [3,4,5] RI- 

; (8) RI ; (8) P 

;[1,2] RI + 

; [3,4,5] RI- 

; (8) RI ; (8) P 

; [1,2] HI + 

[3,4,5] RI- 
(8) RI ; (8) P 

[l,2]RI + 

; [3,4,5] RI- 

; (8) RI ; (8) P 

; [1,2] RI + 

; [3,4,5] RI- 

; (8) RI ; (8) P 

; [1,2]RI + 
; [3,4,5] RI- 



7 1-1 8 1 7 1 

A + A = A e [l,2l 
i-a '-a 1-7 L ' ' 



punch A 

7 1 l ~° 

■ *_ 7 6 [3,4,5] 

7 I 

punch Aj_ 7 and store 

a i-i 7i ei 

A + A = A e [1,2] 

1-7 1-7 1-6 ' 

7 .' 6 1 . . 

- A .-,- A ,-e £ t 3 ' 4 ' 5 ] 



punch A and store for next card 



A i-e + A ,- 6 = A ,. 5 « [1,2] 

7 1 6 1 5 1 

" A H - A ,-e- A .- 5 E t 3 ' 4 ' 5 ! 

5 1 

punch A and store for next card 



4 1-1 5. 1 4 1 

A + A = A e [1,2] 

1-5 1-5 j-4 l ' ' 

7 1 4 1 r i 

-A - ... - A 6 [3,4,5] 

i-7 1-4 

4 1 

punch A and store for next card 

1-4 



3 1-1 4 1 

A + A 

1-4 1-4 



V e [1,2] 

1-3 



7 1 

,_ A 

1-7 



V e [3,4,5] 

1-3 



punch A and store for next card 



2 1-1 3 1 

A + A 

1- 3 I 



2 1 
\ 
1-2 



A . e [1,2] 



7 1 

-A - 

1-7 



V_ 2 € [3,4,5] 



punch A and store for next card 

1-2 



11-1 2 1 

A + A 

1-2 1-2 



7 1 

A 

1-7 



A ,-! « [1,2] 
Al.,6 [3,4,5] 
punch A and store for next card 

store for next card 
[f,+ A , + A + + A ]e [3,4,5] 

L 1 1-1 1-2 -t--..-f j_ 7 J L , , J 
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Note: Wire [3,4,5] exits — [l,2] entries through a Co-Sel picked up at P x . 
Remarks: 



1 i-l 2 i-l 

-f. , - A - A 

i -1 i-2 i- 3 



7 i-l 



A t _ g is stored in [3,4, 5] from previous card. 



i i-i 
A is stored in (1) 

i-2 v ' 



2 i-l 
A i-3 



" (2) 



7 i-l 

A " " " (7) 

i-8 V ' 



Functional Notation: 



i i-i 
A. „ 



2 i-l 

A 

i -3 



2 i 

A 



i-2 



3 A i_l 
A i-4 V 



5 A 1_1 R 

i-6 



4 i-l 



1-4 



3 A' 



'a 1 " 1 

i-a 



8 A i 



V" 1 

i-7 



V 



i-a 



V 



i-7 



i -1 



1 A i 



i-i 



V 



= f 



i-8 i 



[ f + 'a 1 ' 1 

L i-l 



2. i-l 3,i-l 4i-l 5,i-l 6, i-l 7 A i-l -, 

„ + A „ + A +A + A + A +A 

i-2 i-3 i-4 i-5 1-6 i_ 7 i-8 J 



CONTROL RELAY Q 



PILOT SELECTORS 




602-A 




IMMED. P.U. 


SET-UP PULSE 


cl t n o y 


/ RU. 

X OR D D.O 



TRANSFER RELAY 




T N C 



© © 



<Q) POINTS ON PLUGBOARD 
O INTERNAL POINTS 
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NOTES ON 602 -A TIMING CHART 



1) The card feed cycle consists of two machine cycles. On the first of these, a 
card is read at the reading brushes. On the second, the succeeding card is 
read at the control brushes. This second cycle coincides with program 1. 

2) Each machine cycle may be thought of as broken into two sections: 0°thru 120 
and 120° thru 360° . The first section operates control functions and the second 
section consists of the actual calculating operation. Note that all digit times are 
in the second section. 

3) In some cases, entry times shown are not the entire time that the functions will 
receive, but rather the times that they must be impulsed for correct operation. 
Usually, when this condition exists, the relays will pick up but will not hold if 
impulsed at times other than those shown. Also, during multiplication by 2 or 
by 5, the timing of several functions is automatically extended past 312° to 334° . 

4) The "10° " shown under the times for some entry hubs means that the function 
must be impulsed for at least 10° somewhere within the time shown to cause 
proper and safe operation. 

5) Timing of C and CI hubs varies widely during different operations. If it is desired 
to select these hubs, use a selector which is up for 360°. 
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STOP CONTROLS FOR THE IBM TYPE 604 
ELECTRONIC CALCULATING PUNCH 

Bruse Moncreiff 
The Prudential Insurance Company of America 



Stop Controls should be wired so that the error card is the top card in the 
stacker. If one of the controls is incorrectly wired, the operator may look in vain for 
an error in the top card. The machine may have been wired so that the stop occurs with 
the error card still in the machine, to be delivered to the stacker one or two cycles later, 
or with the error card buried in the stacker under either one or two cards. It is important 
that the operator be able to find the error card should an error stop occur. More serious 
than the loss of time is the danger that the operator may decide that since the top card does 
not contain an error, the machine has stopped by accident. 

There are three stop hubs, labeled "1st Reading" , "Punch" , and "2nd Read- 
ing". The entry of an impulse into any of these hubs will cause the machine to stop. The 
difference between them is the number of cycles which are completed after the impulse 
occurs, before the machine stops. If the impulse is wired to "2nd Reading Stop" , the 
machine will stop at the end of the cycle during which the impulse occurs. If the impulse 
is wired to "Punch Stop", the machine will stop at the end of the cycle following the one in 
which the impulse occurs. In other words there is a one cycle delay in stopping. If the 
impulse is wired to "1st Reading Stop", there will be a two cycle delay in stopping the 
machine . 

It should be recalled that five machine cycles are required to move the card from 
the hopper through the machine and into the stacker. The only other facts required to 
correctly wire the stop controls are contained in the 604 Timing Chart. It can be seen 
from this chart that the "Unfinished Program" and "Zero Check" hubs both emit X impulses 
on the cycle following the calculate time in which those conditions are detected. These 
impulses come when the card which caused either of these conditions is at the punch station 
(assuming that the factors used in the calculation were read at "1st Reading"). This impulse 
should be wired to the stop hub which will allow the machine to run until the card in error 
has been delivered to the stacker. An immediate stop would leave the error card just past 
the punches; a one cycle delay would leave the card just past the "2nd Reading" brushes; 
and a two cycle delay would leave the card just delivered to the stacker. The latter is 
the desired situation. Therefore, under the normal circumstances where all factors are 
read at "1st Reading", both "Unfinished Program" and "Zero Check" should be wired to 
"1st Reading Stop". (See chart, line 1, 2. ) 
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If, however, some of the factors are being read at "2nd Reading" and a re- 
calculation and zero check is wired, then "Zero Check" should be wired to "2nd Reading 
Stop". (See chart, line 5.) In cases such as this where factors from two cards (one set 
from "1st Reading" and one set from "2nd Reading") are being operated on during every 
calculate time (604 program) , an "Unfinished Program" signal indicates that both the 
card being calculated and the one being checked must be re-run one cycle at a time. This 
is done by alternately depressing the start and stop keys so that there will be a pause be- 
tween each machine cycle. This gives the electronic unit time to complete the calculation. 
If "Unfinished Program" has been wired to "1st Reading Stop" the two cards to be re-run are 
the top card and the third one down. (See chart, line 8.) If, however, "Unfinished Program" 
was wired to "2nd Reading Stop", the error cards will be the top card and the one to be 
delivered to the stacker two cycles later. The wiring of "Zero Check" may alter this 
situation. It is possible that both lights may come on at the same time. The safest course 
of action is to re-run all the cards in the vicinity and manually check the results. 

Double punched or blank column cards read at either brushes into the "Double 
Punch and Blank Column Entry" hubs are signalled at X time of the following card cycle 
by an impulse from the "DPBC" hub. Double punches or blank columns detected at the 
2nd reading brushes would be signalled in the cycle in which the error card is being deliver- 
ed to the stacker. In order to stop with the error card uppermost in the stacker, "2nd 
Reading Stop" should be wired from "DPBC". (See chart, line 3.) If the card is 
detected at first reading brushes, the "1st Reading Stop" should be wired. 

The product overflow device can be used to signal the presence of a digit 
impulse other than zero from any exit: "1st Reading", "Counter Exit", "General 
Storage Exit", or "2nd Reading". It is most often used to detect a result being read out 
to "Punching" which is greater than the card field assigned to it. The condition arises 
during the cycle in which the error card is being punched, and is signalled by an X impulse 
from "Product Overflow Out" the following cycle. If this impulse were to be wired to 
"2nd Reading Stop", the machine would stop with the error card just past the 2nd brushes. 
The correct wiring would be to "Punch Stop". (See chart, line 4.) This would cause a 
one cycle delay in the stopping, allowing the error card to be delivered to the stacker. 
Similar reasoning can be used to determine the correct wiring when product overflow is 
used for other purposes. (See chart, lines 6, 7.) 

It is important that the operator know what the various error lights indicate on 
each job. Since these controls are pluggable, and their use varies from job to job, the 
meaning of these error lights should be made a part of the operating instructions. These 
instructions should include the type of error (double punch, error in calculation, etc.), 
the location of the error card when the machine stops, the field involved, and the method 
of calculation. This information will give the operator all the information he needs to 
detect and correct errors. This dual attack of correct control panel wiring and full 
information for the operator should prevent errors from slipping through. 
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AN EIGHT-DIGIT GENERAL-PURPOSE CONTROL PANEL 

William P. Heising 
IBM Technical Computing Bureau 



Introduction 

In the course of our work in the IBM Technical Computing Bureau, we have 
found it desirable to set up a 604 control panel to evaluate electronically various 
trigonometric, exponential, and logarithmic series, in addition to the normal arithmetic 
operations, for a CPC or a 604 equipped with punch delay. For these series, a 
number, n , related to the highest power of the argument to be used in the truncated 
series, is read into the 604 with the argument. The series is evaluated by the method 
of continued products on successive program repeats. 

Method 

The eight-digit factors, A and B (with seven decimal digits) are read into 
FS1-2 and FS3-4 respectively. GS3 and GS4 are impulsed to read in on every card 
cycle; a "1" is emitted into GS4 in the high order position, and the number, n, is read into 
GS3 from the card for the series evaluations. The result appears in the electronic 
counter rounded to seven decimal digits (the two lowest positions of the counter are to 
be ignored). The result automatically is read into the electronic storage for factors 
A and B in case this result is to be used as a factor for the succeeding card. 

The operations outlined in Table I require 47 program steps and six calculate 
selectors. The first thirteen program steps are left open to permit additional wiring 
for special operations using the remaining two calculate selectors. These are altered 
to fit the individual problem. Figure 1 shows the programming involved, and Figure 
2 shows the wiring of the calculate selectors. (See pages 80-82) 
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TABLE I 



Operation Code 

or 
Calculate Selectors 




Transferred 


Operation 


4 


A + B 


3,4 


-A + B 


8 


A x B 


3,8 


-A x B 


7,8 


A-J- B 


3,7,8 


-A -r B 


6,7,8 


aTa 


5 


exp (-A) 


3, 5 


exp (+A) 


5,7 


sin A 


5,7 


cos A 


3,5,7 


sinh A 


3,5,7 


cosh A 


5, 8 


log It^aJ 


3,5,8 


log U + A) 


5,7,8 


tanh" A 


3,5,7,8 


tan _1 A 


5,6,7,8 


sin" A 


3,5,6,7,8 


sinh A 



Notes 



B = 


1.2 




B = 


0* 


n 
n 

n odd 
n even 
n odd 
n even 
n 
n 

n even 
n even 
n even 


i 


' 


n even 



*See later section on "Accuracy of Series" 
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Notation 

The planning diagram of Figure 1 uses the first column of program exits for 
"read-out" instructions, and the second column is generally reserved for "read-in" 
instructions. The following notation is used: 



A - FS1-2 = x A xx xxxxx 

- A 

B = FS3-4 = x A xx xxxxx 



Result in electronic counter OOx x A xxxx xxx/xx 



D x = GS3 



n = xxx A 



D 2 = GS4 



10000 
C = GS1-2 ; All assignments are 8-6 



Z.T. 



Ctr.Add 
Ctr. Subtract 
Multiply 
Divide 
Zero Test 



R = Ctr. Reset 

RO = Ctr. Read Out 

RR = Ctr. Read Out and Reset 

MQ = Multiplier Quotient 



Roman Numeral (V) = Suppression Type V 

4 T = Calculate Selector 4 transferred 
15-3 = Third program exit of program step 15. 
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Supp. 


P.S. 


R.O. 


R. I. 




* Notes 


I 


14 


B 


+ 






II 


15 


B 


+ 






II 


16 


RR 


D t 


D 2 RI 




III 


17 


B 


+ 


in 3 rd 


18-2 (- if 3 T) 


IV 


18 


A 


+ * 


* 


18-3 (in 5 th if 7_ T; 


XI 


19 


GS1DO 


GS 2 PU* 




in 3 rd if 7Nand 4T) 


VI 


20 


A 


+ 


GS1PU 


19-2 (unless 7 T) 


VI 


21 


A 








VI 


22 


A 


(GS 2 PU*} 




21-2 and 22-2 (-if 8 T) 


V 


23 


A 


_ * 




23-2 (+ if 3 T) 


VII 


24 


RO 


MQ 






VII 


25 


RR 


C 


out 6 th 




VII 


26 


B 


X 






VII 


27 


C 


MQ 






VII 


28 


RR 


C 


out 6 th 




VII 


29 


C 


+ 






VII 


30 


B 


X 






I 


31 


A 


+ 


in 5 th 




VIII 


32 




1/2 adj. 


in 2 nd 




VI 


33 


RR* 


B 


out 3 rd 


*33-l (RO if 6 T ) 


VI 


34 


GS2 PU 









(continued on next page) 
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Supp. 


P.S. 


R.O. 


R.I. 




IX 


35 


RR 


B 


out 3 


IX 


36 


D 2 


+ 


c th 
in 6 


IX 


37 


GS 2 DO 


GS1PU 




X 


38 


»x 


• 




X 


39 


RR 


C 




X 


40 


C 


+ 


a th 

in 6 


X 


41 


MQ 


c 




X 


42 


D x 


• 
• 




X 


43 


R 






I 


44 


C 


+ 


• o th 

in 6 


I 


45 


MQ 


+ 




I 


46 


B 


- 


in 3 rd 


XI 


47 


D ! 


+ 




XI 


48 




Z.T. 


in 6 th 


XL 


49 


RO 


D ! 




VIII 


50 


Emit 1 


+ 


in 6 th 


xn 


51 


R 


Z.T. 


in 3 


x 


52 


C 


+ 


in 6 


X 


53 


MQ 


+ 




V 


54 


° 2 


+ 


in 6 th 


V 


55 


GS 2 DO 






VIII 


56 


B 


- 


in 3 rd 


None 


57 




1/2 adj. 


in 2 nd 


None 


58 


RO 


B 


out 3 rd 


XIII 


59 


R 


P. Rpt. 




Supp. on 


60 


B 


A 




Non-Zero 











Figure 1. Programming 
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Figure 2. Calculate Selector Wiring 
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Accuracy of Series 

A. Exponentials, Circular and Hyperbolic Sines and Cosines 

The MacLaurin series of these functions converges for all arguments. In 
this general-purpose control panel, the argument and the function must be less than 
ten in magnitude. 

Of the two errors, (a) truncation error of the infinite series and (b) accumu- 
lated rounding error, the latter is the limiting factor on the accuracy obtainable. The 
last term retained in the finite series approximation is Bx n /n! . At n = 40, the 
truncation error is negligible for all values of the argument. The accumulated rounding 
error for arguments of large magnitude, |A| , is proportional to exp |A| . This error 
is 2-4 XlO" 4 for an argument of 9.9. The round-off error for angles near 2 v radians 
is 5xl0~ 6 in the sine and cosine. 

B. Logarithms, Inverse Functions 

These series converge for |A| <1 or |A|« 1 (the latter where the successive 
terms alternate in sign). Because of the slower convergence, a greater number of terms 
is required, in general. The last term retained (or first term omitted if B= 0) is 
Bx n in magnitude (except for arcsin and arcsinh series). For the inverse tangent series, 
if B = JL , the function is quite accurately approximated up to and including arguments 
of unity." Thus for n=990, and B=0. 0005057, the maximum error in the arctangent is 
llxlO -7 (at 44° 58'). For n=90, B=0. 0055688, the maximum error is 1. 7X10" 5 . In 
both these cases, B was chosen so that arctan 1= 7r/4 exactly. In general, the rounding 
errors are negligible, and the truncation error limits the accuracy. The time of 
electronic calculation is approximately 10 program repeats per second. 
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PART 3 



Application 
of the 
IBM Card- Programmed Electronic Calculator 
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INTERPOLATION 

ON THE 

IBM CARD-PROGRAMMED ELECTRONIC CALCULATOR 

Stuart R. Brinkley, Jr. and G. L. Wagner 
United States Department of Interior, Bureau of Mines 



Stirling's interpolation formula for a function y(x), tabulated at equally spaced 
values of the argument x, can be written in the form, 



2 3 4 

y = y o +a t u+a 2 u +a 3 u+ a 4 u + 

where u = (x - x q ) /h, h = x x - x q , 

y = y (x), y = y (x n ), 



(1) 



a, = (1/2) 



A W + A (l) 

-1 o 



-(1/12) 



A (3) + A (3 > 

-2 -1 



(Z) 



(4) 



a 2 = (1/2) A_ x - (1/24) A_ 2 , 



a 3 = (l/12) [C + A^, 



a =(1/24) A 



(4) 



(11 



and where A_ 2 = y (x_ t ) - y (x_ 2 ), 



CD 



A -i = y ( x „> ~ y ( x -i)' 



(2) 



(i) 



id 



(2) 



(1) 

A 



(i) 



- A 
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(3) (2) (2) 

A . = A . - A „ 



(3) (2) (2) 

A . = A - A . 



(4) (3) (3) 

A . = A . - A _ 



In the present application, it is assumed that fifth differences vanish, and that f(x) is 
single-valued. 

Equation 1 is employed for forward interpolation. The table is arranged 
in order of increasing values of the argument. Except at the beginning and end of the table, 
the label subscripts of the table argument are assigned so that 

X « X < X, 

o 1 

If x is less than the third table argument, then the third table argument is taken to be x . 

o 

If x is greater than or equal to the next to last table argument, then the next to the last 
table argument is taken to be x . 

For inverse interpolation, equation 1 is written in the form 
x - x = uh 

o 

u = u o + F(u) (2) 

" = (y - y ) /a t 

2 3 4 

3> 



F(u) = -X (a u +a u + a u ). 



Equation 2 is solved by iteration. The table is arranged in order of increasing values 
of the function. Except at the beginning and end of the table, the label subscripts of the 
table argument are assigned so that 

y (x ) ^ y < y (x x ). 

If y is less than the third table value, then the third table argument is taken to be x . 
If y is greater than or equal to the next to last table value, then the next to last table 
argument is taken to be x t . 
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The circuits for the solution of equations 1 and 2 assume that the quantity y 
is a function of the argument x and of a parameter, a. 

y = y (a, x). 

It is assumed that y is a single -valued, monotone increasing function of x. 

It is assumed that the argument increment h is a one-digit positive number. 

The change in variable 

x ' = a + bx, 

can be employed to transpose into an acceptable form, cases for which y is a single-valued, 
monotone decreasing function of the original argument or for which the argument increment 
h is not a one-digit number. If the function is multiple valued, interpolation except near 
extrema can be carried out on a series of single-valued sections of the table. 

The control circuits are initiated by a control card (first card) identified with 
an X punch. The control card is followed by the table cards in order. No additional control 
cards are employed. Six cards must follow the table card with argument x o . If x o is 
near the end of the table, blank cards in the number necessary to provide six following 
cards are added to the table. For forward interpolation, the control card contains the value 
of the interpolation argument x in the same card field as in the table cards. For inverse 
interpolation, the control card contains the value of the interpolation function y in the same 
field as in the table cards. 

Forward interpolation is performed with Setup Change Switch 2 on. Inverse 
interpolation is performed with Setup Change Switch 2 off. No other wiring changes are 
required. 
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FORWARD INTERPOLATION 
Card Program 





Seq 


Address 


Channel 


Counter 




B C 


B 


c 


2 


3 


4 


5 


6 


1st Card 




00 


1 

y 






(a) 

X 


(a) 

X 






1st Table Card 




00 80 


y.„ 




(a) 

a 


(a) 

-X 






(a) 

y.n 


2nd Table Card 




00 85 


y..« 


y.„ 








(b) 




3rd Table Card 




00 85 


^-n+2 


— n 




J? 




-n 




4th Table Card 




00 85 


y.n+3 


A (d 

-n+l 




A d) 




A<" 

-n+l 




5th Table Card 




00 85 


y.„ + 4 


A (l) 
-n+2 




(d) 
-h 




-n+2 
























nth Table Card 




00 85 


y 


a« 

-2 




-h 




*5 




(n + l)th Table Card 




00 85 


yi 


AW 
-l 




-l dl 








(n + 2)th Table Card 




00 85 


y 2 


^ 




(d) 
-h 




a 1 : 5 




(n + 3)th Table Card 




00 85 


y 3 






JP 








(n + 4)th Table Card 
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y 4 














(n + 5)th Table Card 


2 


00 75 


x-x o -2h 






— 




to) 

-A* 

O 




(n + 6)th Table Card 




00 


y 6 


y-y 








y-y D 




(n + 7)th Table Card 




Punch 






a 




X 


y 





(a) From 3rd Reading 

(b) From Ctr 6 Exit 

(c) From GS 4 Exit 

(d) From Digit Emitter 

Sequence 1 - Differencing Sequence 
Sequence 2 - Interpolation Sequence 
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INVERSE INTERPOLATION 
Card Program 







Address 


Channel 


Counter 






Seq- 


B C 


B 


c 


2 


3 


4 


5 


6 


1st Card 




00 


y 












V 


1st Table Card 




00 85 


y.„ 




(a) 
a 


(a) 
x n 








2nd Table Card 




00 85 


y -„ + i 


y-„-y 








y-y.„ 




3rd Table Card 




00 85 


Y-n+2 


A U1 








- A.„ 




4th Table Card 




00 85 


Y-n+3 


A -n+1 








-n + l 




5th Table Card 




00 85 


y.» + 4 


A W 

-n+2 




( h b) 




- A - n+ 2 






















_ 


nth Table Card 




00 85 


y 


A", 




« 




A° 
-A-2 




(n + l)th Table Card 




00 85 


Yi 


A -1 




( P 




-A-i 




(n + 2)th Table Card 




00 85 


y 2 


A W 

" O 




(b) 
h 




-^ 




(n + 3)th Table Card 




00 85 


y 3 






(b) 
h 








(n+4)th Table Card 




85 


y 4 














(n + 5)th Table Card 


2 


00 73 


y-y -A 










s. 




(n + 6)th Table Card 




00 


y 5 


x-x 




X " X o 








(n + 7)th Table Card 




Punch 






a 


X 






y 



(a) From 3rd Reading 

(b) From Digit Emitter 

Sequence 1 - Differencing Sequence 
Sequence 2 - Interpolation Sequence 
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INTERPOLATION 
604 Program 





Sup 


Sel. 


Factor Storage 


Mult. 
Quot. 


General Storage 


1 2 


I 




No. 


1(6-4) 
3 


2 


4 


Ctr. 


1(6-4) 
3 


2 


4 


F 


1 


1 






RO 


RI 












X 


X 


X 


2 


2 












RO 




RI 




X 


X 


X 


3 


2 










RO 


RI(-) 








X 


X 


X 


4 


3 


IN 




RI 






RO,RC 


















IT 








RI 


RO.RC 








X 


X 


X 


5 


1 
















RO 


RI 


X 


X 


X 


6 


2 










RO 


RI(+) 








X 


X 


X 


7 


2 












RI(+) 




RO 




X 


X 


X 


8 


2 












RO,RC 




RI 




X 


X 


X 


9 


2 


3N 










RI(+) 






RO 












3T 


(EMITTER 2h) 




RI(+)4 








X 




X 


10 


3 




RO 








RI(+) 








X X 


X 


X 


11 


1 












RI(-) 


RO 






K 


X 


X 


12 


1 




RO 










RI 






K 


X 


X 


13 


1 












RO 




RI 




>C 


X 


X 


14 


2 




RI 








RO,RC 








X 


X 


X 


15 


3 


IN 










RI(+) | 






RO 












IT 










RI(-) 






RO J 


t 


X 


X 


16 


1 






RI 






RO 






3 


t 


X 


X 


17 


2 












Ri(+) 






RO 


X 


X 


X 


18 


3 








RO 




RI(-) 






> 


C X 


X 


X 


19 


2 












RO,RC 






RI 


X 


X 


X 


20 


2 




(EMITTER 5) 


RI4 










X 


X 


X 


21 


2 






RO 






MULT(- 








X 


X 


X 


22 


2 






RI 






R05 






1 


x 


X | X 



(continued on next page) 
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Sup 


.Sel. 


Factor Storage 


1 




General Storage 


2 


I 




No. 


1(6-4) 
3 


2 


« 


Mult. 

Quot. Ctr. 


1(6-4) 
3 


2 


4 


F 


23 


2 




(EMITTER 6) 


RI5 










X 


X 


X 


24 


2 








RO 




MULT(+) 








X 


X 


X 


25 


2 








RI 




R05,RC 








X 


X 


X 


26 


2 












MULT(+) 






RO 


X 


X 


X 


27 


2 












RI(-) 5 




RO 




X 


X 


X 


28 


2 












R05,RC 


RI 






X 


X 


X 


29 


2 


3N 
3T 


RO 

(EM] 


[TTER 


h) 




MULT(+) 






RI4 


X 

X 


X 


X 


30 


2 


4N 


RO 








MULT(+) 








X 


X 








4T 


RO 








RI (+) 5 








X 




X 


31 


2 


4N 
4T 








DIV 
DIV 




RO 




RO 


X 
X 


X 


X 


32 


2 












RO, RC 








X 


X 


X 


33 


2 


4N 


RI 






RO 










X 


X 




34 


4 












RI(-) 5 




RO 




X 


X 


X 


35 


4 






RO 






MULT(+) 








X 


X 


X 


36 


4 












R04, RC 






RI 


X 


X 


X 


37 


4 








RO 




RI(-) 6 








X 


X 


X 


38 


4 












MULT(+) 






RO 


X 


X 


X 


39 


4 












R05,RC 






RI 


X 


X 


X 


40 


4 


5T 










RI(-) 6 


RO 






X 




X 


41 


4 












MULT(+) 






RO 


X 


X 


X 


42 


4 


5N 
5T 










R05,RC 
R05,RC 


RI 




RI 


X 
X 


X 


X 


43 


4 


bN 
6N 










MULT(+) 






RO 


X 


X 





(continued on next page) 
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Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


2 


I 




No. 


Sup.Sel. 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


F 


43 


4 


5T 
6T 










MULT(-) 


RO 






X 




X 


44 


4 




(DO GS 2) 




RO 








RI 


X 


X 


X 


45 


5 










DIV 




RO 






X 


X 




46 


4 


6N 










RO.RC 








X 


X 








6T 










R06.RC 


RI 






X 




X 


47 


5 




(PU GS 1) 




RO 


Rl( + ) 








X 


X 




48 


5 




RO 








RI(+)2 








X 


X 




49 


5 










RI 


R02 








X 


X 




50 


5 












RI(-)2 






RO 


X 


X 




51 


4 


6N 


(EMITTER 2) 




RI(-)2 








X 


X 








6T 


(EMITTER 2) 




RI(+) 








X 




X 


52 


5 




(BALANCETEST) 












X 


X 




53 


6 


(PROGRAM REPEAT, PU GS 2) 








X 


X 




54 


4 


7N 
7T 


(EMITTER 4) 
(EMITTER 1) 




RI(+)2 
RI(+)2 








X 
X 


X 


X 


55 


5 




(BALANCETEST) 












X 


X 




56 


7 


( 


PROGRAM REPEAT, PU GS 2) 








X 


X 




57 


4 


7N 










RO,RC 








X 


X 








7T 










RO,RC 






RI 


X 




X 


58 


8 


7N 
7T 


(EMI 


TTER 


h) 




Rl(+) 


RO 




RI 


X 
X 


X 


X 


59 


8 


8N 










MULT(+) 






RO 


X 


X 








8T 








DIV 








RO 


X 




X 


60 


8 


8N 










R02,RC 




RI 




X 


X 








8T 








RO 






RI 




X 




X 
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en 




<u 




s 


z 
o 


Eh 


t— 1 


u< 


H 


CD 


< 




i-l 


c 


o 


d 


Oh 


Si 

Eh 


H 


U 


Eh 


u 


5 






• — l 




0! 




(« 



T3 . 

rs — ' 




1 1 1 1 1 1 1 : 


1 1 - 1 1 1 1 


5" 


CM 

r— 1 


z z z z z z z : 

Eh Eh Eh H Eh Eh Eh : 


Z Z Z Eh Z Z Z 
Eh Z Z Z Z Z Z 


CO 

Si 

o 

■8.2 
0) 


00 

t- 

CO 
CM 


T N N N N T 
N T T N N T 
N N T N N N 
N. N T N N N 
N N T N N N 
N N T N N N 
N N T N N N 


N N T N N N 
N N T N N N 
N N N T N N 
N N N N T N 
N N N N N N 
N N N N N N 
N N N N N N 


CO 

u 

-w o 

o ->-> 

-H O 

CD 

CO 


© 

os 

CO 

t- 

CD 
ITS 

CO 
CM 


TNNNNNNNNNN 
NTNNNNNTNNN 
NNTNNNNTNNN 
NNNTNNNTNNN 
NNNNTNNTNNN 
NNNNNTNTNNN 
NNNNNNTTNNN 


NNNNNNTTNNN 
NNNNNNTTTNN 
NNNNNNNNNTN 
NNNNNNNNNNT 
NNNNTNNNNNN 
NNNNNTNNNNN 
NNNNNNNNNNN 




Si 
d 
U 


1st Card 
1st Table Card 
2nd Table Card 
3rd Table Card 
4th Table Card 
5th Table Card 
6th Table Card 


(n + 2)th Table Card 
(n + 3)th Table Card 
(n + 4)th Table Card 
(n + 5)th Table Card 
(n + 6)th Table Card 
(n + 7)th Table Card 
(n + 8)th Table Card 



o 
a 
u 



T3 

U 
d 

U 

o 

<+1 

Si 
o 

<+1 

w 

CD 

• — i 

o 
>> 

o 



d 

Eh 

CD 
Si 

d 



to 

U 
O 

-t-> 

o 

CD 
CD 



g 



d d 



O 

a 

$1 

CD 



CD 
CO 

Si 

CD 
'f-t 

u 
o 

CO 
CD 
i — i 

u 



CD 
Si 

d 



CO 

u 
o 

-4r> 

o 

CD 
■ — i 
CD 
03 



3 
o 

»— < 

d "-» 

O fi- 
TS ° 
3 CNJ 

,n ■« 
m y 

d cq 

^ CD 
co be 

u a 

O d 

CD ^ 

w 3 

o « 

8s 



3 



d . 

U 35 

■d ° 

3 <M 

^ X! 

d CO 



CD 
bD 
C 

xl 
U 

ft 

3 

CD 
CO 
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INTERPOLATION 
Calculate Selector Circuits 



SEL. 


T 


o 


O 


O R1MQ 


ORI(-) 


O 


1 


N 


-o 


OS1I1 


O RIFS2 


O RK+) 


O 




C 


OSI31 


O SWT 


O PGM 3 


O PGM 15 


o 




T 


O GSI 


O GS2 


O 


OS(-) 


o-s(+) 


2 


N 


O SWT 


O SWT 


O SWT 


OSWT 


OSWT 




C 


-O SB) 


OS{8) 


— 0S(4) 


-*> 


— o 




T 


O EMIT 2h 


ORI4 


OEMIT h 


ORIGS4 


ORI4 


3 


N 


O R0GS4 


O 


OR0FSI3 


O MULT.1+) 


o 




G 


O PGM 9 


O PGM 9 


OPGM 29 


OPGM 29 


O PGM 29 




T 


O R 1 (+) 


O RI5 


O R0GS4 


O 


O 


4 


N 


O MULT(+) 


O 


O R0GSI3 


OR0MQ 


O RIFSI3 




C 


O PGM 30 


O PGM 30 


OPGM 31 


O PGM 33 


O PGM 33 




T 


ORI(-) 


O RI6 


O R0GSI3 


ORIGSI.3 


OMULT(-) 


5 


N 


o 


O 


O 


ORIGS4 


OMULT(+) 




C 


O PGM 40 


O PGM40 


OPGM 40 


O PGM 42 


O PGM 43 




T 


O R0GSI,3 


O RIGSI3 


ORIH-) 


O 


OSWT 


6 


N 


O R0GS4 


O 


oRI(-) 


ORI2 






C 


O PGM 43 


OPGM 46 


O PGM 51 


O PGM 51 


OS(7) 


7 


T 
N 


O EMIT 1 
O EMIT 4 


O RIGS4 
O 


O Rl(-H 
O EMIT h 


O R0GSI3 


OSWT 


ORIGS4 


1 




C 


O PGM 54 


OPGM 57 


OPGM 58 


O PGM 58 


OS(6) 


8 


T 
N 


O DIV 

O MULT{+) 


O ROMQ 
O RO.RC 


O 


O 


O SWT 


O R02 


O 


1 




C 


O PGM 59 


OPGM 60 


OPGM 60 


O 


O S(5) 



Note: PGM means program step. 
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53 


3 


O 




3 


U 


O 


-M 


ft 


o 


tf 


p— ( 


w 

Eh 


CD 



5 3 



cu 



CD 



Q 

z a 

CM CC 




= "!! 

CL « 



o <o 

o _) 

3 UJ 

a. en 
o 



o * 
a. en 



-, id 
5; <n 




o eg 
°_i 
z> "J 





_J 
cr> 



« N 
3 UJ 

a o 
o 



o 




en 




c 



111 u u 
« in . 
O I z 




<l - 

O J 



2 




rt 



e> i- 
a. in 





2 
o 



95 



INTERPOLATION 
Co Selector Circuits 



SEL 


T 


^j o PISEL n CTR 6<+) 
3-2N 


CTR5(+) 
CTR SRO 


Q PU CASEL 
3-8 


1 


N 


PI SEL 
°5-2 N ° ° 















C 


oCTR3(+) OCTR3H O— i 




OPISEL 


OCCY 






T 




— — o o 




OCTR 6(+) 


OCTR 4(+) 




2 


N 


O 


— o o 




O 


O 






C 


OCCY 


OACY O 






O 






T 


O 


— o o— 




OCTR 2(+) 


O 


3 


N 


O— 




— — o o 


O 


O 




C 

T 


o 




OPUCO OCCY 


OCCY 


O 








PISEL 
3-2 C 


o 








4 


N 


o 


o o 




O 


o 




C 


o 




p> 






OCCY 


o 




T 


o— 




o(5) 


0<7) 




0(3) 


OCTR 3 NBT EX 


5 . 


N 


— o 


o(3) 


O(B) 




0(5) 


OCTR5NBT EX 




C 


OC AOD(I) 


s* 


1 




? 


OPISEL 7-1 C 




T 


o 


y 


o 


0(8) 


I " " 


O 


6 


N 


o 




PU FLD 
SEL 


O(0) 


o(0) 


O 




C 


o 




o(D 


OB ADD(2) 


OB ADD(I) 


O 


7 


T 
N 
C 






o 


0(7) 


O 
O 

o 


CTR 5(+) 
R0GS4 






— O O OCADDI2) 


— o 




T 
N 
C 








8 


o oth) o o 


o 











Note: CCY means card cycle 
ACY means add cycle 
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9 
8 
7 
6 
5 
4 
3 
2 



INTERPOLATION 
Field Selector Circuit 



TABLE FUNCTION 3rd READING 
o O O o o o o o 




C 



CHANNEL B 
o o o o o o o 
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INTERPOLATION 
Digit Circuits 



2nd Reading: Parameter a - Comparing Entry (1) 

1st Card X - PU Pi Sel 1 
3rd Reading: Argument x - CTR 4 (5-1) Entry 

- Co Sel 8T (5-1) 
Parameter a — Comparing Entry (2) 

- CTR 2 Entry 
Function y - Field Sel C (6-1) 

Field Sel 1 (6-1) - Channel B (6-1) 
Channel B <6-l) - FS1, 3 Entry 

- CTR 6 (6-1) Entry 

- CTR 5 (6-1) Entry 
Digit Emitter (h) - Co Sel 8N (4) 

Co Sel 8C (5-1) - CTR 3 (5-1) Entry 
GS2 (5-1) - Channel C (5-1) 
GS4 (5-1) - CTR 5 (5-1) Entry 



98 



INTERPOLATION 
Differencing Sequence 





Sup. 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No. 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


R 




y 2 


(2) 


A< 2) 
A -2 


A (3) 
A -3 


A (3) 
A -2 


\ 


a'." 


*? 


1 


1 




RO 


RI 

A (2) 












4 


3 








RI 

A (3) 
-2 


RO,RC 








5 


1 














RO 


RI 


10 


3 


RO 








RI(+) 

y 2 








11 


1 










RI(-) 

A w 

i 


RO 






12 


1 


RO 










RI 

y 2 






13 


1 










RO 




RI 

A<? 




15 


3 










RI(-) 

A 2 ) 








RO 


16 


1 




A (2) 

A o 






RO 








18 


3 






RO 




RI(-) 

A (3, 








P 




y 2 


A <2) 
o 


*-"! 


A (3 * 

"-2 


*Z 


y 2 




A 

O 
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FORWARD INTERPOLATION 
Interpolation Sequence 





Sup 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No. 


1(6-4) 
1 3 


2 


4 


1(6-4) 
3 


2 


4 


R 




x-x o -2h 


a!' 


(2) 

A-i 


A (3) 
A -2 


13) 


y 2 


A W 





2 


2 










RO 




"? 




3 


2 








RO 


RI(-) 
24a 4 








4 


3 




RI 
24 a 4 






RO,RC 








6 


2 








RO 


RI(+) 

A ( i 2 








7 


2 










Ri(+) 
12 a 3 




RO 




8 


2 










RO,RC 




RI 
12a 3 




9 


2 


(EMIT! 


T2R 2h) 






RI(+) 4 
2h 








10 


3 


RO 








RK+) 

x-x 








14 


2 


RI 

x - x <, 








RO,RC 








15 


3 










Rl(+) 

A u1 








RO 


17 


2 










M(+) 

2A? 






RO 


18 


3 






RO 




RK-) 
2(a 1+ a 3 ) 








19 


2 










RO,RC 






RI 
2(a t + a 3 ) 


20 


2 


(emit: 


:er 5) 




RI4 
5 










21 


2 




RO 






MULT(-) 
-12a 4 









(continued on next page) 
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Sup. 


Factor Storage 


1 




General Storage 


No. J 


1(6-4) 
3 


2 


4 


Mult. 
Quot. 


Ctr. 


1(6-4) 
3 


2 


4 


22 


2 




RI 

-12a 4 






R05 








23 


2 


(EMIT 


rER6) 




RI5 
6 










24 


2 






RO 




MULT(+) 
12 a 2 








25 


2 






RI 

12 a 2 




R05,RC 








26 


2 










MULT(+) 
12(a 1+ a 3 ) 






RO 


27 


2 










RI(-)5 
12 ax 




RO 




28 


2 










R05,RC 


RI 

12 a 1 






29 


2 


(EMIT 


TERh) 












RI4 
h 


30 


2 


RO 








RI(+)5 

X " X o 








31 


2 








DIV 

u 








RO 


32 


2 










RO,RC 








34 


4 










RI(-)5 
-12 a 3 




RO 




35 


4 




RO 






MULT(+) 
-12(a 3 + 
ua4) 








36 


4 










R04,RC 






RI 

-12(a 3 + u a 4 ) 


37 


4 






RO 




RI(-)6 
-12 a 2 








38 


4 










MULT(+) 
-12(a 2+a3 u 

2 &\ 

+ a 4 u ) 






RO 



(continued on next page) 
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Sup 


Factor Storage 


4 


Mult. 
Quot. 


Ctr. 


General Storage 


No. 


1(6-4) 
3 


2 


1(6-4) 
3 


2 


4 


39 


4 










R05, RC 






RI 

-12(a 2 + a 3 u 
+ a 4 u^ 


40 


4 










RI(-)6 

-12 a l 


RO 






41 


4 










MULT (+) 
-i2(y u -3fcV 






RO 


42 


4 










R05,RC 


RI 

-12(y-y„) 








u 




43 


4 










MULT(-) 

12(y-yJ 


RO 






46 


4 










R06,RC 


RI 
I2(y-y ) 






51 


4 


(EMITTER 2) 






RI(+) 
2 








54 


4 


(EMIT 


TER1) 






RI(+) 2 
12 








57 


4 










RO,RC 






RI 
12 


58 


8 










Ri(+) 
12(y-y c ) 


RO 






59 


8 








DIV 

y-y„ 








RO 


60 


8 








RO 






RI 

y-y 




P 




x-x 


-12 a 4 


12 a 2 


y-y D 




i2(y-y ) 


y-y„ 


12 
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INVERSE INTERPOLATION 
Interpolation Sequence ) 





1 
Sup. 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No. 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


R 




y-y -a? 


A? 


A.? 


A-? 


A-? 


y 2 


A? 


^ 


2 


2 
GS1 










RO 




A -i 




3 


2 

GS1 








RO 


Rl(-) 
24 a 4 








4 


3 

GS1 




RI 

24 a 4 






RO,RC 








6 


2 
GS1 








RO 


m(+) A <* 

A 2 








7 


2 

GS1 










RI(+) 
12 a 3 




RO 




8 


2 
GS1 










RO,RC 




RI 
12 a, 

3 




9 


2 
GS1 










ra(4 V 

-A 






RO 


10 


3 

GS1 


RO 








RI(+) 

y-y„ 








14 


2 
GS1 


RI 

y-y 








RO,RC 








15 


3 

GS1 










A o 






RO 


17 


2 

GS1 










m(+) a> 
2^„ 






RO 


18 


3 

GS1 






RO 




RI(-) 
2(a l + aa) 








19 


2 

GS1 










RO,RC 






RI 
2(a 1+ a 3 ) 


20 


2 

GS1 


(E 


MITTER 


5) 


RI4 
5 










21 


2 
GS1 




RO 






Mult (-) 
-12 a 4 









(continued on next page) 
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Sup 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


22 


2 

GS1 




RI 

-12 a 4 






R05 








23 


2 

GS1 


(EMIT 


'TER 6) 




RI5 
6 










24 


2 
GS1 






RO 




MULT(+) 
12 a 2 








25 


2 
GS1 






RI 

12 a 2 




R05.RC 








26 


2 
GS1 










MULT(+) 
12(a 1+ a 3 ) 






RO 


27 


2 
GS1 










RI(-) 5 
12 a t 




RO 




28 


2 
GS1 










R05,RC 


RI 
12 a x 






29 


2 
GS1 


RO 








MULT(+) 

6(y - y ) 








30 


2 
GS1 


RO 








MULT(+) 

I2(y-yj 








31 


2 

GS1 








DIV 
u 




RO 






32 


2 

GS1 










RO,RC 








33 


2 

GS1 


RI 






RO 










34 


4 










RI(-)5 
-12 a 3 




RO 




3b 


4 




RO 






MULT(+) 
-12(a 3 +a 4 u) 








36 


4 










R04,RC 






RI 
-12(a 3 +a 4 


u) 


37 


4 






RO 




RI(-)6 
-12 a 2 











(continued on next page) 
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Sup, 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No, 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


38 


4 










MULT(+) 
12 ai F/u 2 






RO 


39 


4 










R05,RC 






RI 

12 a x F/u 2 


41 


4 










MULT(+) 
12ajF/u 






RO 


42 


4 










R05.RC 






RI 

12 a x F/u 


43 


4 










MULT(+) 
12 a t F 






RO 


44 


4 


(DOGE 


»2) 




RO 








RI 
u 


45 


5 








DIV 
F 




RO 






46 


4 










RO,RC 








47 


5 


(PUGS 


1) 




RO 


RI(+) 

F 








48 


5 


RO 








RI(+) 2 
u ' 








49 


5 








RI 

u' 


R02 








50 


5 










RI(-) 2 
u'-u 






RO 


51 


4 


(EMITTER 2) 
1 






RI(-)2 
u'-u-2(10" 4 ) 








52 


5 


(BALANCE TEST) 












53 


6 ( 
NEG 
3AL 


PROGRAM 


1 REPEA 


T, PUGS 


2) 










54 


4 


(EMITTER 4) 






RI(+) 2 
u'-u+-2(10" 4 









(continued on next page) 
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Sup. 


Factor Storage 


Mult. 
Quot. 


Ctr. 


General Storage 


No. 


1(6-4) 
3 


2 


4 


1(6-4) 
3 


2 


4 


55 


5 


(BALANCE TEST) 












56 


7 
POS 
BAL 


(PROGRA 


MREPE 


AT, PUG. 


3 2) 










57 


4 










RO,RC 








58 


8 
GS2 


(EMIT r 


fERh) 












RI 
h 


59 


8 
GS2 










MULT(+) 

X " X o 






RO 


60 


8 
GS2 










R02,RC 




RI 

X " X o 




P 




U o 


-12 a 4 


12 a 2 


u 


,„ 


12 a t 


X ' X o 


h 
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